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ABSTRACT. We consider a classical system of point particles interacting by means of 
a short range potential. We prove that, in the low-density (Boltzmann-Grad) limit, the 
system behaves, for short times, as predicted by the associated Boltzmann equation. This 
is a revisitation and an extension of the thesis of King [9] (appeared after the well known 
result of Lanford [10] for hard spheres) and of a recent paper by Gallagher et al [5]. Our 
analysis applies to any stable and smooth potential. In the case of repulsive potentials 
(with no attractive parts), we estimate explicitly the rate of convergence. 
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1 Introduction 



In a well known paper in 1975, O. Lanford presented the first mathematical proof of the 
validity of the Boltzmann equation for a system of hard spheres, for a sufficiently small 
time. The starting point was the series expansion describing the time evolution of the 
statistical states of a hard-sphere system. This series is the solution of a hierarchy of 
equations formally established by C. Cercignani in 1972 [2], following previous ideas due 
to H. Grad [6]. 

The main idea of Lanford is to compare such a series expansion with the one arising 
from the solution of the Boltzmann equation, claiming the term by term convergence in 
the so called Boltzmann-Grad limit (BG limit in the sequel). The restriction to short 
times is due to the fact that the two series have been proven to converge absolutely only 
for a small time interval. Actually it was remarked in [20] that the Lanford's approach is 
a Cauchy-Kowalevski kind of argument. 

In [10], although all the main ideas, as well as the strategy of the proof, were clearly 
discussed, the details were missing. The complete proof was presented later on in [9], [16], 
[19], [18] and [3]. 

We mention also that the ideas of Lanford can be apphed to derive the Boltzmann 
equation globally in time, in the special case of an expanding cloud of a rare gas in the 
vacuum [7, 8]. 

Shortly after the appearance of the Lanford's paper, F. King in his unpubhshed thesis 
[9] approached the same validity problem for a particle system interacting by means of a 
positive, smooth and short range potential. In this case the basic starting point was not 
the usual BBGKY hierarchy, but a variant of that due to H. Grad [6] (we shall call it the 
"Grad hierarchy" in the sequel) making the system more similar to a hard-sphere one. 
More precisely, in [6] only the first equation of this hierarchy was discussed, while the full 
hierarchy was introduced and derived in [9]. 

The Boltzmann equation considered by King was written in unusual form. Namely, 
calling / = f{x,v,t) the distribution function, 

{dt + v-Vx)f{x,v,t)= / dvi du {v-vi)-u\f{x,v[,t)f{x,v',t)-f{x,vi,t)f{x,v,t)\ 

(1.1) 

where S*^ = {u G 5*^1 {v — vi) ■ u > 0}, S'^ is the unit sphere in M^, (^,^1) is a pair of 
velocities in incoming collision configuration -see also [1]- and {v', v[) is the corresponding 
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pair of outgoing velocities defined by 

{v' — V — uj\uj • (v — Vi)] 
^ ^ . (1.2) 

v'l = Vi +uj[uj ■ {v - Vi)] 

Here u = oj{v, V) is the unit vector bisecting the angle between the incoming relative 
velocity V = vi—v and the outgoing relative velocity V = v[ — v' as specified in the figure 
below. 
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Figure 1: The two-body scattering. We denote by p the impact parameter expressed in 
microscopic unities (p G [—1,1]) and by x = x(P) 1^1) f^e scattering angle (% G (— 7r,7r], 
X > in the figure), while © is the angle given by the relation x — ~ 20. We call 
scattering vector the function uu — uj{i', V). 

A more handable and usual form for the Boltzmann equation is obtained by expressing 
everything in terms of oj, namely 

{dt + v-V,)f{x,v,t)^ [ dvi ! du B{uj,V){f{x,v[,t)f{x,v',t)- f{x,vi,t)f{x,v,t)} 

(1.3) 

where B{uj,V)/\V\ is the differential cross-section of the potential under consideration 
with respect to the solid angle u (in the case of hard spheres, the two formulations (1.1) 
and (1.3) coincide since ou — u). 

After many years, the argument has been recently reconsidered by I. Gallagher et al 
in a long and self-contained paper [5] pointing out some important facts, surprisingly not 
discussed in the previous literature. In particular, the term by term convergence is not 
innocent because B is, in general, not bounded and even not defined as a single-valued 
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function. For instance, for smooth positive and bounded potentials (considered by King 
himself), u ^ u is not globally invertible and B is unbounded. Therefore the difficulty is 
that one has to exclude concentration of measure on certain small sets in the phase space 
leading to an evolution much different from the typical Boltzmann behavior. These "bad" 
events are: (i) the long time two-body scattering; (ii) the recollisions, i.e. the presence of 
a given pair of particles undergoing two or more collisions. The latter is the main obstacle 
in proving that the particle system behaves as predicted by the Boltzmann equation. 

In [5] the authors prove the validity of the Boltzmann equation under the hypotheses 
that the potential is well behaving in this sense, namely that the cross-section exists as a 
single-valued and sufficiently regular function. In the present paper we show that, under 
very general assumptions on the potential, the Boltzmann equation can indeed be derived 
in the form (1.1). We review the results in [9], completing some parts of the proof and 
taking care of some inconsistencies. Once the Boltzmann equation has been derived in the 
form (1.1), the passage to the form (1.3) is a matter of analysis of the two-body problem. 
If the cross-section is not a single-valued function, the function B appearing in (1.3) can 
be still expressed as a sum of the contributions arising from each monotonicity branch. 

The approach discussed in [5] makes explicit use of the cross-section as a tool for the 
control of recollisions. In the present paper the aim is to establish a proof that docs not 
depend on any detail of the scattering process. In particular, the term by term convergence 
(which is the most delicate point in the proof of our main results) is treated in a different 
way from the one in [5] and [9]: see Section 7.1 for a presentation of the problem, and 
Sections 7.2, 7 for a quick abstract and an exphcit constructive proof respectively. In 
our method a very useful tool is a tree expansion describing the time evolution of the 
marginals of a statistical state. This is presented in Section 6. In Section 2 we introduce 
the mechanical system of particles under examination and make some preliminary remark 
about it, while in Section 3, along the lines of [9], we derive the Grad hierarchy, that is 
the starting point of our study. In Section 4 we fix the hypotheses on the initial data and 
state our main results. In Section 5 we present the uniform short time estimates on the 
series expansion for the evolution of the marginals. The results in Sections 3 and 5 are well 
known by [9] and [5]: we discuss them here briefly for the sake of completeness. Finally, 
in the Appendix we give sufficient conditions on the interaction for having a bounded or a 
single-valued cross-section. 

One advantage of the methods developed in this paper is that they allow an explicit 
estimate of the error in the convergence to the Boltzmann equation, as soon as one has 
explicit estimates of the interaction time of the two-body process in the space of the 
scattering parameters. Moreover, convergence is established in a strong sense, that is 
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uniformly outside a precise pathological null-measure subset of the phase space. 

The analysis of sections 2-7 can be applied to any smooth and repulsive potential, 
enlarging the class of interactions considered in [5]. If the potential has an attractive part, 
there is also an additional difficulty due to long time scattering phenomena and to the 
presence of trapping orbits in the two-body process. For the sake of clearness we treat this 
case separately in Section 8, where we explain how the proof can be adapted to extend the 
convergence result, assuming the stability of the interaction. 



2 The hamiltonian system 

We consider a system of N identical classical point particles of unit mass, moving in the 
whole space and interacting by means of a two-body, short range potential We denote 
by (gi, Vi, - ■ ■ , qiy, vjy) a state of the system, where Qi and Vi indicate the position and the 
velocity of the particle i, and Qiir) the position of particle i at time r. The A?"— particle 
Hamiltonian is 

JV N 

^=2^^'+2E^(^^-^^)- (2-1) 

1=1 i,j=l 

The dynamical flow is obtained by solving the Newton equations 

= (2-2) 

where F{qi — qj) — Fij — — V$(gi — qj) is the force due to the particle j, acting on the 
particle i. We will assume $ to be smooth enough in order to have existence and uniqueness 
of the solution to (2.2) for any initial datum such that qi ^ qj (see Hypothesis 1 in Section 
2.2, and Section 8). 

Consider now a small parameter e indicating the ratio between the macro and the micro 
unities. We pass to macroscopic variables defining 

X = eq; t = st . (2.3) 

In these variables the equations of motion become 

^^^\t) ^^^f( ^'^^) ~^^'(^) ^ (2 4) 

In order to have a kinetic picture, a tagged particle, say particle 1, must deliver a 
finite number of collisions in a macroscopic unit time. As a consequence, the density Ne^ 
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must vanish. More precisely N should be 0{e~'^). Indeed, assuming the characteristic 
interaction length of the potential $ to be one in microscopic variables, namely $(g) = 
if \q\ > 1, consider the tube spanned by particle 1 in the (macro) time 1: 

{^\^^y^-^iit)\<e] . (2.5) 

The number of particles in the tube is the number of particles potentially interacting with 
particle 1 in a macroscopic unit time. Hence, if N — 0{e~'^), such a number is expected 
to be finite. Therefore the scaling we will consider is 

N ^oo, e^Q, Ne^ = l'^ > 0, (2.6) 

for a system of N particles obeying to (2.4), where I > will be proportional to the mean 
free path and will be fixed to one for notational simplicity. 

The scaling (2.6) is usually called low-density hmit and it is equivalent to the BG hmit 
originally introduced for the hard-sphere system, [6]. If we want to picture the dynamics 
in macroscopic variables, we can say that a triple collision - namely a situation in which 
three or more particles are simultaneously interacting - will be very unlikely. Moreover a 
two-body collision - namely a scattering process involving only two particles - will take 
place on a scale of time of 0{e), but since the force is 0{e~^) it will produce a finite effect. 
In other words the expected dynamics is qualitatively similar to that of the hard-sphere 
systems. 



2.1 Statistical description 

We want to describe our system from a statistical viewpoint. We will use bold letters for 
vectors of variables, for instance 

will be the notation for the state of particles 1, ■ ■ ■ ,j and j + 1, ■ ■ ■ , j + n respectively, 
having position and velocity {xi,Vi). As usual we introduce the phase space 



Mn = i zat e 



t,6N 



\xi-Xk\>0, i,k^l---N, kj^i^. (2.1.2) 



Consider a probability distribution with density W'^ , which is initially (and hence 
at any positive time) symmetric in the exchange of the particles. Its time evolution is 
described by the Liouville equation, which reads as 

{dt + CN)^ ^0, (2.1.3) 
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where the Liouville operator jCn is 



with: 



(2.1.4) 



i=l 
1 ^ 

= -Y.^^^r^v, (2.1.5) 



and F,,,- = -V$ (^) . 

Remark. For simpUcity wc shall assume in this subsection, as well as in Section 3 below, 
that is a smooth (say C^) function of its variables over all A^at x M"*", with Vi ■ VxiW^ , 
-fij • V^.Vr^ e L^{Mn)- The assumption is used to write the Liouville equation as a partial 
differential equation and to perform partial integrations. This is not really needed to state 
our results: we will weaken the regularity hypothesis later on by using a density argument 
(see Proposition 1 and the discussion before it). 

We introduce the marginals gij^(zj,t) of the time evolved measure {zN,t), defined 

by 

gf{zj,t) = J dzj,j,_jW''{zj,zj,N_j,t) , (2.1.6) 

which denote the probability distributions of the first j particles (or of any other fixed 
group of j particles). Clearly — . 

Prom (2.1.3) and (2.1.6) it follows that the family {gj^}jLi satisfies the well known 
BBGKY hierarchy ([6]): 

(d^ + ±v.. gf ^\±F (^) . V..,f 

- E / dx,^, J dv,^,F [^.1^^^ . . (2.1.7) 

Notice that, for a fixed j, the interaction term in the left hand side of Eq. (2.1.7) is, in a 
sense, negligible because the collisions among a tagged group of particles are unlikely (the 
potential is indeed vanishing as soon as e is smaller than \xi — Xk\)- Moreover the integral in 
the right hand side is 0{e^). The right hand side, which is due to the interaction between 
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the group of the first j particles with the rest of the system, is whenever N — 
which is exactly the reason why we perform the low-density scaling. However, instead of 
using the above hierarchy, not very well suited for such a scaling, we will introduce, in 
Section 3, another set of equations. 



2.2 The two— body scattering 

Let us discuss here the scattering process between two particles, which will play a crucial 
role in what follows. We turn back to microscopic unities, where the potential is assumed 
to have range one. 

Let qi,vi,q2,V2 be positions and velocities of two particles which are performing a 
collision. It is well known that this two-body problem can be reduced to a central-motion 
problem setting the origin in the center of mass: 

— ^ — = 0, q^qi-q2- (2.2.1) 

Then the evolution is given by 

0(r)=2F(g(r)) . (2.2.2) 

The above equation of motion is "almost" explicitly solvable. In particular, fixed the relative 
velocity V = V1—V2 (hence fixed a value of the energy in the center of mass), one can restrict 
his attention to the control of the scattering function u = u^u). Since the scattering takes 
place in a plane, this amounts to control the function = 0(p) (see Fig. 1 in Section 1). 
The classical integral formula expressing in terms of the modulus of the incoming relative 
velocity \V\, the potential $ and the impact parameter p will be written in the Appendix 
(see Eq. (A.2)). That formula is not so easy to handle with, so it will not be employed in 
the present section. 

In what follows it will be rather crucial to have an estimate on the scattering time r*, 
namely the measure of the time interval for which \q{T)\ < 1. To this purpose, we need to 
state our precise assumptions on the potential. 

Hypothesis 1 The two-body potential $ = $(g),g G M^, is radial, with support \q\ < 1 
and not increasing in \q\. We assume either $ G C^(M^), or $ G C^(M^ \ {0}) and $(5) 
+00 as q ^ 0. 

The smoothness assumption is needed to ensure existence and uniqueness of the flow evo- 
lution for the system of N particles, while the monotonicity is introduced to allow a simple 
control on the scattering time r*. We defer more general cases to Section 8. 
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We will use sometimes the notational inconsistency for which $(r) = ^||g|=r-- 
Consider the central motion given by Eq. (2.2.2) with the initial conditions describing 

the two particles just before the interaction, namely g(0) S'^, q{0) — V and \V\ > 0, 

V -ly <0. Denote 

L^\u AV\^\pV\ e[0,V] (2.2.3) 

the magnitude of angular momentum, being p the impact parameter (Fig. 1). A rather 
general estimate on is the following: 

Lemma 1 Under Hypothesis 1 it is 

n<- (2.2.4) 

for some ^4 > 0. 

Proof. From the conservation laws one derives the well known formula expressing r* as 
a function of V and L : 

/•I I 

n = V2 dr— ^ , (2.2.5) 

where r* is the minimum distance from the origin, r* = inf^g(o,T,) !?('?") I) related to V and 
L by 

T = ^ + 2$(r.). (2.2.6) 

The effective potential, i.e. the potential of the reduced one-dimensional motion (which 
is the evolution of the radial coordinate in the system of the center of mass), is the 
L— dependent function 

2^e//(r) = 1^ + 2$(r) - ^ , r e [0, 1] . (2.2.7) 



We can write 



1 

dr 



- ($e//(r*)-$e//(r))^/' 

1 1 

< r / dr . . (2.2.8) 

(min[o,i](-$;,,))^-/- 

Denote improperly $' the derivative with respect to r of the function Since $' < 
and 

Kffir) = $'(r) - 1^ , (2.2.9) 
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the result follows easily. ■ 

The estimate in Lemma 1 tells us that = 0{{pV)^^). This has the advantage to 
be general and sufficient to our purposes. Clearly the bound can be improved in many 
cases. Singularities in the scattering occur whenever the colhsion is central (p = 0) and the 
energy corresponds exactly to a point of vanishing force {V^/A = $(r),$'(r) = 0). This 
kind of singularities does not exist if the potential is unbounded at the origin and strictly 
repulsive: for instance for potentials diverging at the origin with a power law, formula 
(2.2.4) can be easily replaced by r* = 0{V~^). Prom (2.2.8) it can be noticed also that the 
singularity for low energies may appear only if $ goes to zero smoothly (C^) in r = 1. 

We conclude by introducing a map which encodes all the properties of the two-body 
interaction. The scattering operator I is defined over 



where co — 00(1/, V) is the scattering vector, see Fig. 1. It follows that V ■ v — —V'-v'. In 
particular, V • u' > 0, i.e. I sends incoming to outgoing configurations. 

The following property of I will be used in the validation of the Boltzmann equation. 

Lemma 2 I is an invertihle transformation that preserves Lehesgue measure. 

Proof. Of course, being the dynamics reversible, , V) — V) (see Fig. 1) and 
is defined in the same way as X. 

To see that X is measure preserving, we fix cartesian coordinates on the plane where 
the scattering occurs, and call (j) the angle formed by V and the first axis (with (j) growing 
when V rotates counterclockwise), a the angle formed by V and v and such that sin a is 
the impact parameter p (with the convention a G [7r/2, (37r)/2], see Fig. 1). Restricting to 
the plane of the scattering, we have the parametrization V = {\V\ cos0, \ V\ sin0), u = a. 
In the variables 0, a the action of X is simply described by: 




(2.2.10) 



by: 



X{u, V) = (^', V) 




(2.2.11) 



V = {\V' \ cos (j)',\V'\ sine/)') 



v — a 
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where 

\v'\ = \v\ 

< a' ^n-a . (2.2.12) 

0' = </> - x(sinQ;, |\^|) 

v 

Note that a' G [— 7r/2, 7r/2]. The first equation is conservation of energy, the second conser- 
vation of angular momentum, and the third holds by definition of scattering angle (Fig.l). 
It can be shown that x is a differentiable function of its arguments (see the discussion in 
the Appendix). Moreover, the determinant of the jacobian of the transformation (2.2.12) 
has modulus one, independently of the form of %. This concludes the proof. ■ 



3 The Grad hierarchy 

In this section we derive a hierarchy of equations for a family of quantities which are 
very close to the marginals introduced in the previous Section 2.1. This allows to put the 
dynamical problem in a form somehow similar to the one arising in considering hard-sphere 
systems and more suitable for the study of the low-density limit. 
We define the reduced marginals 

f^{zj,t) = / dz,. ^_,-iy^(z„ z,, t) , (3.1) 

where 

5(xj) = |z = (x,'t;) e >£ for all A; = l,--- ,j| . (3.2) 

It is clear that the functions /j^, for any j, are asymptotically equivalent (uniformly on 
compact sets in Aij in the BG limit) to the usual marginals. 
Consider now a configuration z^r = (zj,Zj jv-j) such that 

\xe - Xk\> e (3.3) 

for all £ = 1, ■ ■ ■ , j and k = j + 1, ■ ■ ■ , N. Since the range of the interaction is e, the 
interaction between the group of the first j particles and the rest of the system is vanishing. 
Therefore the Liouville equation (2.1.3) on such configurations becomes: 

dtW + CLW"" + ll]W'' + Ll^W"" = , (3.4) 

where is defined as in (2.1.5) with the sums running from j + 1 to N . 
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As already said in the Remark at page 6, we make at the moment the regularity 
assumptions needed to justify Eq. (3.4) and all the steps in the following derivation (see 
the mentioned Remark). 

Integrating Eq. (3.4) with respect to dzj^N-j over S{x.j)^~^ we obtain: 



f dzj,^_jv, ■ V,,M/^(z^, t) , (3.5) 

i=i Js{^jr-i 

where we used that 

/ civ,,^,_,4^W^^ = . (3.6) 
The first sum is handled by the divergence theorem yielding 

- / dzj+i ■■■ da(xi)dvi ■■■ dz^ivi ■ Vi)W^ , (3.7) 

, 1 J SUA JdSU-i) J SUA 



i=j+l 



where is the outward normal to S{x.j) and da{xi)dvi is the surface measure. Note that 
if Zk e dS{xj), there exists an index k e {!,••■ :j} such that \xk — Xi\ — e. Moreover, 
if Zj e A4j, there is only one such an index for almost all Xi with respect to the surface 
measure da{xi). Hence dS{xj), as regards the x-dependence, is the disjoint union of pieces 
of spherical surfaces. We call such pieces crjk(xj), that is 

j 

dSi^ij) = y (7jfc(xj) X (3.8) 

k=l 

where 

(Tfe(xj) C {x I \x — Xi\ — e} . (3.9) 

We set = |'^'°~^'| . Using the symmetry of we have N — j identical integrals, for 
which (3.7) becomes: 

-(^-i)y] / da{xj+i) dvj+i{vj+i ■ iyk,j+i) dzj+i^N-j-iW^ . (3.10) 

The integration domain of the last integral in (3.10) is not S'(xj_|_i)^^''^-'^, as it would 
be necessary to recover fj^^ and close the equation. We could reduce this integration to 
S{'Xj^i)^~^~^ , and this would produce a small error in the BG limit. Nevertheless, we 
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want to establish an exact and closed equation, for which we need a further work. Let 
i — {ii, • ■ ■ , i^} be a subset of ordered indices of {j + 2, • • • , N}, with ii < ^2 < • • • < im- 
We put Zi — {zi^, • ■ ■ , Zi^}, and introduce the set 

Ai(x,_|_i) := < Zi C S(xj)"^ such that, for each £=!,••• , m, min \xk — Xi,\ < e> . 

I feeiu{j+i} * J 

(3.11) 

A generic configuration in S^Xj)'^"^"^ differs from S'(xj_|_i)^~^~^ because some particle, 
say particle ii, could overlap with particle j + 1, this meaning that \xi^ — Xj+i\ < e. If 
this is the case, we consider the maximal cluster of overlapping particles with indices i. 
Definition (3.11) gives the subset of such cluster-configurations. The other particles are 
far apart the group with indices 1, • • • ,j,j + 1,*, therefore each of them is in S{x.j+i,Xi). 
Then, we can decompose the integration domain S{xj)^~^~^ in (3.10) as a union of disjoint 
sets to obtain 



j N-j-l 



k=l m=0 i:|i|=„-'<^fe(xi) 

(3.12) 

JAi(xj+i) 

where we denoted |i| the cardinality of i and we used (3.1), as well as the symmetry of 
W^, to integrate out the not-clusterized variables. All the terms in the are identical 
so that the result can be written 



-E E iN-j)iN-j-l)---iN-j-m) I da{x^+,) 

fe=l m=0 -^f^fcCxj) 

dVj+l{Vj+i-iyk,j+l) [ ^^^^^/ili+m(zj+i+m) , (3.13) 



/ 



where 

Am{'Xj^i) := ^Zj^i^rn C S'(xj)™' such that, for each £ = j + 2, ■ ■ ■ , j + 1 + m, 



mm 

,j+l+m 



^ \xi — Xj(\ < . (3.14) 
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The second sum in Eq. (3.5) is 

j p j 

+{^-j)y2 da{xj+i) dvj+i{vi ■ Uij+i) dzj+i^N-j-iW^{zN,t) . 

i=l J<^i{^j) -'kS Js(x,)^-^-i 

(3.15) 

Repeating the same step we did before to go from (3.10) to (3.13), we readily arrive to the 
following hierarchy of equations (which we call Grad hierarchy): 

N-j-l 

{dt + Cj) /f = E AVi+^/ili+m, 1<J<N, (3.16) 

m=0 

where the operator Cj = depends also on s through its interacting part (2.1.5), and 
A^J+^+mf!^l+mi^J, t) = {N- j){N -j-l). ..(N-j-m) 

•VsM |^,+^s_^,|>e}(z/) / dVj+i{Vj+i - Vi) ■ 

/j^i+^(zj-,a;i + i/£,'yj+i,Zj+i,^t) , (3.17) 



/ 

Ja, 



with Xj+i — Xi + i/e in the argument of A^. We indicate with !{.}(•) the characteristic 
function of the set defined by the condition in the curly brackets. 
In particular it is 

■{Vj+i - V,i) ■ Z//j^i(Zj, Xi + PS, Vj+i,t) , 

= £2(7V-j)q+i/i:i(z„t), (3.18) 

where (3.18) defines Cj_^_i, which is the same collision operator appearing in the hard-sphere 
case, see [10]. Actually it is clear that, in the BG limit, this term is the only 0(1) term in 
the sum in the right hand side of (3.16). Indeed, for m > and fixed j, the size of A^j_^l_^.^ 
will be 

0(7V'"+^£V") , (3.19) 

the £^"* coming from the successive integrations in the domain A^(xj+i). This means that 
we are in a situation quite similar to that of the hard-sphere system [10], and we can hope 
to derive the Boltzmann equation in a similar manner. 
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3.1 Series solution 



Consider the dynamical flow obtained by solving the Newton equations (2.4) for a system 
of j particles: 

^W.lg^(flW^), (3.1.1) 

where i and k run now from 1 to j. Denote by Tj{t)zj the solution of this system of equations 
with initial datum Zj. The action of this flow on the functions is given by the interacting 
flow operator Sj{t), defined as 

S^(t)g(zj) = g(r.(-t)zj) . (3.1.2) 

We may represent the solution of Eq (3.16) by means of a perturbative expansion, that 
is just the iteration of the Duhamel formula: 

//W = E E dti dt2--- din (3.1.3) 

n=0 mi,...,m„>0: 

■Sl{t - il)^^+i+^j<S|+i+„^(il - ^2) • ■ ■ Aj^^^Y,7=imi^j+n+j:7=imii'tn)ff+n+j:'^^^mii^) ' 

where fj^{t) = /j^(-,t), and /j^(0) are the reduced marginals of the initial probability 
distribution. This expansion will be our main tool. 

We derived Eq. (3.1.3) assuming sufficient smoothness of the initial distribution (see 
Remark on page 6). However, by using a density approximation, (3.1.3) can be proven to 
hold for a general class of initial measures. The argument can be found in [19] page 281, 
or [17] page 18 for cases of hard-sphere dynamics, and it can be applied also to general 
smooth potentials. We list the main steps in what follows. A different approach based on 
a weak formulation may be found in [5] . 

Consider the collection of integration variables in the right hand side of (3.1.3), which 
we call A = {ti, ■ ■ ■ jinji^i, - • • , i^n, Vj+i, • • • , fj+n; ^ij+n+^.mj) 866 also (3.17). The reduced 
marginal in the integrand takes a form /j^^^^n ^^.(yj+n+x;r-i"»i(^j' ■^)' 0); with Zj e A4j 
and A in the integration domain (for the understanding of the detailed structure of the 
map, (zj. A) — >■ yj^n+^7_j^mi, we defer to the discussion in Section 6 of this paper). Now, 
let us write the expansion in (3.1.3) for a generic measurable probability density W^. To 
have a nice control on the integration over large velocities, we shall assume the exponential 
decrease < c^e~^^i=^^i for some c, ^ > 0. Since y is a Borel map (as follows directly 
from measurability of the partial mappings {zj,t) — )■ T|(t)zj), the expansion makes sense 
for the reduced marginals of , and the integrals therein are absolutely convergent. 
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To recover identity (3.1.3), we use that there exists a sequence of smooth densities 
W'^''^ which evolve according to (3.1.3), satisfy the exponential bound, and approximate 
: limW^^'''' = a.e. on M.n- Since the densities evolve according to W^{z]^,t) — 

7— )-0 

W^(TU-t)zN) (and same equation for W^^^), we also have \imW^^^(t) = W^(t) a.e. on 
A^TV and, consequently, lim/^''^(t) = f^{t) a.e. on Aij. We are left with the problem of 
taking the limit of the right hand side of (3.1.3). Using the measure preserving property 
of the flows T^(i), it can be shown that y is not singular, in the sense that y^-*^ maps null 
sets in j_|_„_|_j^^_^ to null sets of values of {zj, A) in its domain. This fact, together with 
the gaussian estimate, allows to apply dominated convergence, thus concluding the proof. 
Summarizing, we have the following 

Proposition 1 Let W'^ be the density of a probability distribution on A4n with reduced 
marginals . Suppose that is Borel measurable, symmetric in the particle labels and 
such that < c^e~^^i^i'"i for some c, /3 > 0. Then the evolved measure at time t > has 
reduced marginals f^{t) given by Eq. (3.1.3), for almost all points in A4j. 

Remark. It is important to observe that deflnitions of the operators A'j and C| (re- 
spectively (3.17) and (3.18)) involve a trace problem, so that they arc well posed if they 
act over functions which arc at least continuous over a.a. points of the spheres of center 
Xi and radius e (see deflnition (3.17)). Nevertheless, this is not relevant to our purposes, 
since we will work only with operators of the form J dsAjSj{s). These last are indeed 
well deflned over functions f^ satisfying the hypotheses of Proposition 1, by virtue of the 
nonsingularity of the map . 

For future convenience, let us conclude this subsection by giving some more definition. 
The subscries associated to the dominant term of (3.1.3) (that with all = 0) defines a 
new sequence of functions which we call ■ 

fj{t)^y2al{j) dtj dt2--- dtn (3.1.4) 

n=0 

■5J(t - ti)Cj_^i5J+i(ti - ^2) ■ ■ ■ Cj+„>SJ^_„(tn)//tn(0) ■> 

al{3) := £2"(7V - j)[N - j - I) ■ ■ ■ {N - j - n + I) , (3.1.5) 

where we used definition (3.18). Notice that in the BG limit Q;^(j) = 0{1). 

Finally, it will be convenient to decompose the collision operator C?^]^ in the following 
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form: 



k=l 



- Vj+i) • i/|/j^i(zj-,a;jk + i/£,Vj+i,i) (3.1.6) 



where 



S\^{v\ {vk - Vj+i) • > 0} , 

Sl^{iy\ {vk - vj+,) • 1/ < 0} . (3.1.7) 
3.2 The Boltzmann hierarchy 

In this subsection we treat formally the solution to the Boltzmann equation (1.1) as we did 
in Section 3.1 for the interacting system of particles and compare heuristically the results. 
Let / be a solution to Eq. (1.1). Consider the products 

fjizj, t) = fitf^izj) = fiz,, t)fiz2, t)... fizj, t) . (3.2.1) 

It is easy to show that the fj solve the hierarchy of equations 

{dt + £;) = C,+i/,+i, 1 < j < 00 , (3.2.2) 

where 

j 

Cj+i = Ck,j+i (3.2.3) 

k=l 

Ck,j+1 = ^k,j+l ~ ^k,j+l 

t)^ dudvj+i{vk - Vj+i) ■ vfj+i{zi, ■■■ ,Xk,v'^,-- - , Zj,Xk, t^^+i, t) 

^kj+ifj+ii^j^ t)^ diydvj+i{vk - Vj+i) ■ iyfj+i{zi, ■■■ ,Xk,Vk,--- , Zj,Xk, Vj+i, t) 
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and 

Vj+l = Vl + UJ[UJ ■ {Vk - Vj+i)] 

CO = C(;(z/, — Vk) being the scattering vector (see Fig. 1). 

The infinite hierarchy of equations (3.2.2) (which does not express nothing else than the 
Boltzmann equation whenever the factorization property (3.2.1) holds) is called the Boltz- 
mann hierarchy. Proceeding as before, we may represent its solution by the perturbative 
expansion around the free flow: 

rt rti rtn-i 

fM = dt2--- dtn 

n>0 

■Sj{t — ti)Cj+iSj-^.i{ti — ^2) • • ' Cj+nSj+n{tn)fj+n{0) , (3.2.5) 

where now Sj{t) is the free flow operator, defined as 

Sj{t)g{zj) = g{xi - vit, vi,--- ,Xj- Vjt, Vj) . (3.2.6) 

Note that: 

- Eq. (3.1.3) is an identity which expresses (well defined by means of the iV-particle 
flow) in terms of a finite sum of operators acting on the initial sequence //^(O); 

- fj^, Eq. (3.1.4), is just a technical definition; 

- Eq. (3.2.5) is a series whose convergence must be proven. 

As for the hard-sphere case [10], it is possible to show that such a series is indeed convergent 
for a short time. We will show it in Section 5. 

Remark. This last result implies also local existence and uniqueness of the solution to 
the time-integrated version of the Boltzmann hierarchy in the class of continuous functions 
such that fj{t) < c'e~'^^»=i for some c, (3 > (see e.g. [3]). In particular, in the case 
of initial product states, factorization is propagated in time and each factor is the local 
solution of the time-integrated Boltzmann equation. 

Reminding the discussion at the end of Section 3 and the fact that jCj equals zero for 
s small, we shall guess that (3.2.2) is what one gets just letting e go to zero in the Grad 
hierarchy (3.16). To do so, assume for simplicity that /j^^ is continuous. Then, we may 
try to rewrite the action of C^.'^^^ on f^-^ in such a way that the function is evaluated in 
incoming collision states. Call = £t* the time of interaction of particles k and j + 1. 
Since the scattering process is, in macroscopic variables, almost instantaneous, we assume 
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that the other particles do not interact in the same time interval. By the continuity of the 
flow it will be 




(3.2.7) 



hence 




52 xM3 



diydvj+i\{vk - Vj+i) ■ u\fj^^{zi, ■■■ , 



^k-i Vf^, • • • , Zj, Xk, t) 



(3.2.8) 




I dudvj+i{vk - Vj+i) ■ iyf^+i{zi, ■■■ , 

5ix]R3 



, Xk, Vf,, ■ ■ ■ , Zj, Xk, ) 



where in the second step we simply changed u ^ —u. 

We stress that this heuristic discussion is somehow dangerous. In fact, the required 
continuity property of /j^i, even when true for any fixed N, is lost in the limit. This is 
why we work with integral formulas instead of partial differential equations. The rigorous 
version of the above (standard) argument, which will be presented in Section 7, resorts 
to the convergence of (3.2.5) to (3.1.3), and requires only continuity of the limiting initial 
data fj{0). 

4 Assumptions and results 

We establish here the hypotheses under which we will work. Beyond Hypothesis 1 on the 
interaction potential stated in Section 2.2, we assume 

Hypothesis 2 The initial condition for the Boltzmann equation is a continuous probability 

density fo over M.^, satisfying the bound 



for some /3 > 0. 

Moreover, indicating by H{zj) the j— particle Hamiltonian written in macroscopic vari- 
ables. 



sup 6 fo{x,v) < +00 



(4.1) 




(4.2) 



we have 
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Hypothesis 3 The initial probability density on A4n is a Borel function Wq symmetric 
in the particle labels. Its reduced marginals /oj, J = ,N satisfy the bound 

fU^y^^^^ ^ (4-3) 

for some a > 0. 

By Proposition 1, the /(^- are good initial data for the evolutions (3.1.3) and (3.1.4). 
Note also that Hypothesis 3 implies that we are fixing correlations even at time zero. 
Indeed, if the interaction potential diverges at the origin, f^j{zj) — > exponentially when- 
ever Xk ^ Xi ior k ^ i. Therefore, initial product states are excluded. This situation is 
similar to that of hard-sphere systems, in which an overlapping of any pair of particles 
is not allowed. Even if the potential is bounded, but positive at the origin (which is the 
case of stable interactions), product states are forbidden by Hypothesis 3. In fact, near 
the diagonal (xk = Xi) the factor e^^^^^^ can grow exponentially with j^. 

Our last hypothesis is 

Hypothesis 4 The family f^j converges to /oj := f^'' as N ^ oo, uniformly on compact 
sets in M.j. 

We are now ready to state our first result. Let /j^(i) be the reduced marginals at time 
t, evolved according to Eq. (3.1.3) and let fj{t) be defined as in (3.2.1) and (3.2.5) (which 
will be proven to be an absolutely convergent series in Section 5). Define also the subset 
of particles that cannot collide pointwise under the free evolution: 

Qj = |zj e Mj s.t. {xi - Xk) A {vi - Vk) 0^ . (4.4) 

Theorem 1 Under the Hypotheses 1-4, there exists to > such that, for < t < to and 

hm ff{t) = f,it) (4.5) 

uniformly on compact sets in Qj. 

Theorem 1 is formulated and proven in the same spirit of [10] and [9]. As we shall see 
in Section 7.2, the proof, based on geometrical arguments, is abstract and does not give 
informations on the rate of convergence. However, the result can be improved under 
quantitative assumptions on the rate of convergence and the continuity of the initial data, 
as explained in what follows. 
Define 

Mj{5) = eR^^ \xi-Xk\> 5, i,k = l---j, ki^i^ (4.6) 
for (5 > 0. We assume 
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Hypothesis 5 For some C > 0, 

sup e§ |/3_^.(z^.) - /„^^.(z,) I < {Cye . (4.7) 

Moreover, for some L > 0, 

ei^' \fo{x, v) - fo{x', v)\<L\x-x'\. (4.8) 
Then we have the following: 

Theorem 2 Assume the Hypotheses 1 — 5. Then, for all Zj E f2j,j > and t sufficiently 
small, there exists a positive Eq = £o{zj) such that, for s < Sq and Ns^ = 1 

\ffizj,t)-fj{zj,t)\<C^e\ (4.9) 

where C > 0, 7 > are suitable constants. 

Observe that Hypotheses 4, 5 are a natural notion of convergence compatible with 
continuity of /o and estimate (4.3) (which prevents convergence on the diagonals Xi — Xk). 
To clarify this point, we construct some explicit example in the next subsection. 

4.1 An example of initial datum 

In the following we present a sequence of reduced marginals satisfying Hypotheses 3-5. Set 

W^o^(zAr) = ^/o^^(z^) n l{|..-x,|>e}(Ziv) , (4.1.1) 

l<i<k<N 



= / dz^f^^'iz^) n 1{|..-..I>.}W (4.1.2) 



where 

"^^^^ l<i<k<N 

is the partition function and /o is some density satisfying Hypothesis 2 and (4.8). The 
reduced marginals are 

^ l<i<k<j 

with 

F^'izj) = [ dz,,^/o^(^-^)(z,,^) (fl fl l{|.,_.,|>,}(z^) 

■ I n M\-^--.\>e}M] . (4.1.4) 
\j+l<i<k<N / 
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Let us estimate F^{zj)Zj^^. First observe that, for some Co > 0, 

ZN_j{l - CoNey <Zn< Zn-j ■ (4.1.5) 

The upper bound is obvious consequence of the normahzation of /q. As regards the lower 
bound, note that 

Zn = dzN-lfo''^~^\zN_i) TT l{\xi-Xk\>s}{'^N-l) 

J^'^"-'^ l<^<k<N-l 
» N-1 

■ / dZNfo(zN) TT l{\xi-XN\>s}(^N) 
= / ciZjv-l/o^^^"^^(ZiV-l) TT l{\xi-Xk\>e}{zN-l) 

• / dZNfoizN) 1 - ^{\xi-xj,\<e}izN) 

>Z^_,{l-Co{N-l)e^) , (4.1.6) 

for instance taking Co = (47r/3)||/o||oo- Eq. (4.1.5) follows by iteration. We can also show 
that 

The upper bound is immediate, while the lower bound follows from 

F^(z,) > / dz,,^/o^(^-^)(z,,^) ( 1 - E E l{l^.-.l<4(^iv) ) 

n l{|xi-xfc|>£}(Zj,Ar) j 
<i<k<N J 



.j+l<i< 

3 . 



> ZN-j-j{N -j)Coe'ZN-j-i , 

Zn-j 



>Zr,-Al-3NCoe'^^] , (4.U 



noticing that (4.1.5) implies ZN_j_iZ~\. < {I-CqNs^)-^. If Ne^ = 1 and N is sufficiently 
large. Equations (4.1.5) and (4.1.7) give in turn the bounds 

1 - 2Coj£ < < ^ (4.1.9) 

Zn (1 - Cqs)-' 



and, in particular, 

(4.1.10) 



Zn N-^oo 



21 



uniformly in Zj £ Aij. 

Now it is easy to check that, for N sufficiently large, the Hypotheses 3-5 are verified. 
Hypothesis 3 follows from Hypothesis 2. Hypothesis 5 (hence 4) follows from the estimates 
in (4.1.9). 

In definition (4.1.1) for the initial density we could also replace the product of charac- 
teristic functions by e~'^^*<'=*( see [5]. This defines a sequence of states which are, 
in a sense, the maximally uncorrelated states for which the Hypotheses are satisfied. 

Finally, other families of initial data exhibiting a slower rate of convergence (and im- 
plying possibly a slower convergence in Theorem 2) can be easily constructed, for instance 
enlarging the cut-off in (4.1.1). If in formula (4.1.1) e is replaced by e"'' with 7' e (2/3, 1], 
then, proceeding as before, we obtain 

sup |/o,(z,) - /i;.(z,)| < (C")%-^+^^' . (4.1.11) 
4.2 General strategy of the proof 

The proof follows the main ideas of [10], adapted to the present context. The vahdity 
argument is based on a comparison among the series for the A?^— particle system (3.1.3), 
and the Boltzmann series (3.2.5). 

- First, we prove that both the expansions are absolutely convergent series, for suffi- 
ciently short times and uniformly in e : see Section 5. As a consequence of the estimates in 
Section 5, it follows also that (3.1.3) and (3.1.4) are asymptotically equivalent in the BG 
limit. 

- Then, it remains to prove the term by term convergence of (3.1.4) to (3.2.5). To do 
this, it is necessary a preliminary detailed analysis of the generic term. This is presented in 
Sections 6 and 6.1 for the series (3.1.4), and in Section 6.2 for the Boltzmann series (3.2.5). 
The structure of the generic term is described with the help of a convenient representation 
of formulas in terms of tree graphs. It turns out that any given term can be expressed as 
an integral over a set of special backwards-in-time trajectories of clusters of particles. 

- The proof of the term by term convergence is carried out in Section 7, using in a crucial 
way the picture introduced in Section 6. The issues arising from the convergence will be 
first discussed in Section 7.1, while the abstract proof leading to Theorem 1 and the explicit 
estimates leading to Theorem 2 will be presented in Sections 7.2 and 7.3 respectively. 
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5 Short time estimates 



The aim of this section is to prove that, for times t smaller than a certain to, the expansion 
for fj^{t), Eq. (3.1.3), can be bounded uniformly in e. The Boltzmann series solution Eq. 
(3.2.5) turns out to be an absolutely convergent series for the same values of t. Moreover, 
the difference between fj^{t) and fj^{t) (defined by (3.1.4)) is negligible in the limit. These 
results arc easily established by assuming the bounds on the initial data in Hypothesis 2 
and 3. Here we will follow [9] (sec also [5, 19, 18, 3]). 

To begin with, we notice that our assumptions on the initial data make natural the 
introduction of the norms: 

11^,11^ = supe'^^^^^^l^.^z,)! , gj-.Mj^R, /3 > , 

||f/||/3,a = supe"°^||c/j||^ , c, = {^.}~^, a>0. (5.1) 

By the energy conservation 

\\S'{t)g\\f,,a = \\9h,a , (5.2) 
for all /3 and a for which the right hand side makes sense. 
The crucial estimate is the following: 

Lemma 3 Let gf : Aij — >■ R &e a sequence of continuous^ functions with g^ — Q forj > N 

Sm>0 ^i+l+mfi'i+l+m f 
_ - J j=l 

Then, for (5' < (5 and a' > a there exists a pure constant C > such that, for e small 
enough, 

^ ^ [vw^m^) ^ ■ 

Proof. From definition (3.17) we find 



IS2 

2 "j+l 



r dz 



^Thc continuity here is required only for simplicity of notation, since it assures well posedness of the 
operator action: see the Remark after Proposition 1. If that is not true, the lemma must be reformulated 
for ds f^^_^y A^S^{s) (that is what we really need to control for the proof of Proposition 2 below). This 
can be done in an obvious way using Eq. (5.2) and adding a factor t"/n! in the right hand side of the 
estimate. 
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Here we used that e'^{N — j) < 1 and the positivity of the interaction (Hypothesis 1), for 
which 

_^ j+l+m 

i/(z,+i+„) = H{z,) + i/(z,,i+^) > H{z,) + - ^l- (5-5) 
The last integral in the right hand side is bounded by 

(5.6) 



/3 y V 3 

so that (5.4) is smaller than 

h'^^.^mUCf^eTY. \ ^^^+1(1^^1 + k.+i|)e-'^^^-'''e-i^.V , (5.7) 

where we used again positivity of the interaction, and is a suitable constant. The 
Cauchy-Schwarz inequality implies 



X^d^^ii + if^j+ii) < 



i=l 



3 

\ i=i 

which inserted into (5.7) leads to 

Summing over m and taking the supremum over j with weight e~"'^ we readily get (5.3) 
for £ smaller than a constant depending only on /9, a. ■ 

Let us apply Lemma 3, together with (5.2), to the right hand side of (3.1.3). We 
proceed by iteration. For a given n > 0, we partition the intervals [/3/2,/3] and [q;,2q;] in 
n intervals of the same length ^ and ^, and then apply the above results n times. The 



outcome is 



n>0-^0 ^0 Jo 

+n 

n>0 

<\\fo\kaY.^tC'p,<^T ^ (5-10) 



n>0 



for suitable constants Cs q , C'^^^, having used Stirling formula in the last step. Hence 
we obtained a geometric series which converges for t sufficiently small (and the radius 
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of convergence is explicitly computable in terms of the other constants). Now the same 
argument can be applied in a straightforward way to the series (3.1.4) and (3.2.5). Thus, 
we have proven the first statement of: 

Proposition 2 In the Hypotheses 2 and 3, we have absolute convergence of the series 
(3.1.3), (3.1.4) (uniformly in the BG limit for e small enough) and (3.2.5), for allt < to — 
to{P,a). Moreover, for some C" > 0, if e is small enough, 

lirW-/^WII/3A2a<C"e. (5.11) 

Proof We just need to prove Eq. (5.11). Set C^g^ = \ e^{N - j)q^-^gf+-^ \ . With 
the notations of Lemma 3 and proceeding in the same way, we observe that 

ii(^^ - ng^y,.' < supe--'^- uu^ragf+.^jy 

■'— m>l 

< supe-'^- Y.W'\\3f^^^^h +^-) 

<C''e{^ ^ +^^\\\g''\\Ba (5.12) 

for suitable C^'^ > 0, having used (5.9) in the first inequality, and e sufficiently small in 
the second. Therefore, proceeding as in (5.10), 

writ) - r{t)h/2,2a < E ^ E (%'ic,,annmw 

n>0 ■ k=l ^ ^ 

n>0 

which gives the result with C" depending on (3, a and on the initial datum. ■ 



6 The tree expansion 

In the proof of Theorems 1 and 2 it is convenient to represent each term of the expansions 
(3.1.4) and (3.2.5) as more explicit integrals of the initial data, f^j and /oj respectively. 
As we will see in the present section, it is natural to express such terms by means of binary 
trees which help us to visualize the various contributions. 

Consider first Eq. (3.1.4) which, reminding Eq. (3.1.6), we rewrite as 



^^w = E<(^')EE*(-i)i'^"| / dtJ'dt,... f 



dtr 



n=0 (T„ k„ 

■Sj{t - tl)C^;'|"j+i5|+i(ii - ^2) • • • Cl'^j+n^j+n{in)fo,j+n : (6-1) 
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where 

n 

Cn = (c^l, • • • ^Crn), CTj = ± , knl^^CTjl, 

i=l 

k„ fel=l fc2=l fcra = l 

We introduce the n— collision, j— particle tree graph, denoted r{j,n), as the collection 
of integers ki, - ■ ■ ,kn that are present in the sum (6.2), i.e. 

ki e Ij, k2 e Ij+i, ■■■ ,kn& Ij+n , with = {1, 2, • • • , s}, (6.3) 

so that we shall write 

E* = E ■ (6-4) 

kn r(j,n) 

Note that the number of terms in the sum is j{j + 1) • • • (j + n — 1). The name tree graph is 
justified by the fact that it has a natural graphical representation. This is best explained 
by an example: see Figure 2 which corresponds to r(2, 5) given by 1, 2, 1, 3, 2. 




Figure 2: Tree graph r(2, 5) = 1, 2, 1, 3, 2. We have also drawn a time arrow in order to 
associate times to the nodes of the trees: at time ti the hne j + i is "created". Lines 1 and 
2 exist for all times; they are called "root lines". 

Given a tree graph r(j, n), and fixed a value of <t„ and of all the integration variables 
in the expansion (6.2) (times, unit vectors, velocities), we can associate to it a special 
(e— dependent) trajectory of particles, which we call interacting backwards flow (IBF in 
the following), since it will be naturally defined by going back in time. The rules for 
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the construction of this evolution will be explained in the next section. The notation for 
a configuration of particles in the IBF will make use of Greek alphabet i.e. C^{s), where 
s e [0, t] is the time and there is no label specifying the number of particles. If s e (i^+i, ^r) 
(with the convention to — t, tn+i — 0) we have j + r particles: 

= (Cl(s),-- - ,C1+.(S)) e Mj+r for S e {tr+l,tr) , (6.5) 

with 

Q(s)-(m,vKs)), (6.6) 

the positions and velocities of all the particles being respectively 

The reason to introduce these trajectories is that we want a more explicit expression 
of each term of the expansion (6.1), namely our purpose is to write Eq. (6.1) as 

N-j 

where 

/n 
dA{K, u^, v„„) n B'{'^^■, vj+i - vUU))fo,j+n{Cm , (6-9) 
i=l 

dA is the measure on M" x 5*^" x M^" given by 

dA(t„, Vj^n) = ^{ti>t2 ->tr,}dti . . . dtndui . . . dv^dvj+x . . . dvj+n , (6.10) 
and we use the short notation 

B%Ui;vj+i - r]l.{ti)) = \vi ■ {vj+i - 77fe.(ii))|l{a,^..(^^._^,_^e^(t.))>o}l{|^|_^.(ti)-^|(ti)|>£ vfc^fcj ■ 

(6.11) 

In other words, in the generic term 7^^ (zj, t), the initial datum f^j^^ integrated, with the 
suitable weight, over all the possible time-zero states of the IBF associated to r(j, n), cr„. 

6.1 The interacting backwards flow (IBF) 

Let us construct C^{s) for a fixed collection of variables r(j, n), an, zj, t„, u^, ^j,n, with 

t = to>ti>t2> ■■■ >tn> tn+l = , (6.1.1) 
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and Un satisfying a further constraint that will be specified below. The j root lines of the 
tree graph are associated to the first j particles, with states Ci; ' " " Xj- Each branch j + i 
{£ — 1, • • • ,n) represents a new particle with the same label, and state Cj+i- This new 
particle appears, going backwards in time, at time in a collision state with a previous 
particle (branch) ke e {1, ■ • ■ j — with either incoming or outgoing velocity according 
to a£ — — or — + respectively. 

More precisely, in the time interval {tr,tr-i) particles 1, • ■ ■ ,j + r — l flow according 
to the usual dynamics T^_,_^_i. This defines Cj+r-ii-'^) starting from Cj+r-ii'^r-i)- At time tr 
the particle j + r is "createcf' by particle kr in the position 

C^+ritr) = + UrS (6.1.2) 

and with velocity Vj+r- This defines C^(tr) = (Ci(^r),-- - Xj+ri'^r))- After that, the evo- 
lution in (tr+i,tr) is contructed applying to this configuration the dynamics T^_,_^ (with 
negative times). The characteristic function in the collision operator (3.1.6) (or the second 
characteristic function in (6.11)), is a constraint on Vj. implying that no third particle is 
closer than e to the pair kr,j + r at the time tr- 

We have two cases. If (Tr = — , then it must be Ur ■ (fj+r — Vlri^f)) < 0. In this case the 
velocities are incoming and no scattering occurs, namely after tr the pair of particles moves 
backwards freely with velocities i]{tr) and fj+r- If o"r = +, we require i^r' {vj+r ~Vkri^r)) ^ 
so that the pair is post-collisional. Then the presence of the interaction in the flow T^_,_^ 
forces the pair to perform a (backwards) scattering. The two situations are illustrated in 
Fig. 3. 

Remark. It is very important to note that between two creation times tr,tr+i any pair 
of particles among the j + r, different from the couple {kr,j + r), can possibly interact 
by reaching (or having from the beginning) a distance smaller than e. These interactions 
arc called recollisions, because they generally involve particles that have already interacted 
at some creation time (in the future) with another particle of the IBF. In our language, 
recollisions arc the "interactions different from creations". Though recollisions are expected 
to be unlikely, we will have to analyze them with special care, since they are the main 
responsible of the different behavior of the particle dynamics from the Boltzmann evolution. 

6.2 The Boltzmann backwards flow (BBF) 

The discussion of the two previous sections can be repeated, with minor changes, for the 
case of Boltzmann series (3.2.5). The interacting backwards flow is now substituted by the 
Boltzmann backwards flow (BBF) <^(s) , for which we use the same notations of (6.5)-(6.7) 
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Figure 3: At time tr, particle j + r is created by particle kr, either in incoming (cTj. = — ) 
or in outgoing (cr^ = +) collision configuration. Particle kr is called progenitor of particle 
j + r. 

with the superscript e omitted. The BBF is introduced exactly as the IBF, see Section 
6.1, except for the following differences: 

- the interacting dynamics is replaced by the simple free dynamics; 

- in the right hand side of (6.1.2) the second term is missing, i.e. the created particle 
appears at the same position of its progenitor; 

- there is no constraint on Vj. other than the one implied by the value of cr^; 

- a ar — +, to determine the state of particles in {tr+i,tr), before applying free evolution 
we have to change velocities according to {r]kr{tr),Vj+r) — >■ {r)k^{t~),rjj+r{'tr))j where — >■ 
denotes the scattering rule depicted in (1.2) and Figure 1. Here rikX^t) indicates the hmit 
from the future, and r]k^{t~) the limit from the past. 

Eq. (3.2.5) can then be rewritten: 

oo 

M^j^t) = E E ' (6.2.1) 

n=0 r{j,n) CTn 

where 

/n 
dA(t„, v,-„) H B(uf, vj+i - rjkXtt))fo,j+n(m) , (6.2.2) 
i=i 
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and 



Remark. An important point concerning the Boltzmann backwards flow deflned in the 
previous section is that, in a given interval {tr+i,tr), the velocities r]i{s) = r]i{t~), besides 
being constant, depend only on the velocities of particles in the future of the BBF, and on 
the vectors of impact i^i, • • • , i^r, but not on the interaction times ti, - • • ,tr- This simple 
structure of the BBF will be exploited later on (see for instance Equation (7.2.11)). 

In the proof of the term by term convergence it will be used a change of variables 
transforming integrals over outgoing variables into integrals over incoming variables. This 
is simply the scattering operator of (2.2.11) applied to an interaction of the BBF. We 
conclude this subsection introducing such operation. 

Fix r(j, n), 1 < r < n, Vj e R^-^ and deflne the transformation X^^^ = '^v-'r{jn) '■ 

I^"^^ {^n, Vj,n) = {l^r-l, ^r,n-r, ^j,r-l, V^, Vj+r,n-r) (6.2.4) 

where only the r— th couple (ur^Vj-^-r) is changed according to 

ly'r = -'^r + "^UJriuJr ■ l^r) for Z/^ • {Vj+r - Vkritt)) > 

K = l^r for ■ {Vj+r - Vkr{tt)) ^ " (6.2.5) 

K = Vj+r{tr) - Vkritr) 

Here Ur = oji^r, Vj+r - rjkri'tt))- 

Lemma 4 The transformation X^**^ is a one-to-one, measure preserving map. 
Proof. X^^) is the composition of the two transformations: 

{Ur,Vj+r) > {l^r, K) 

Vr = Vj+r - Vkrit^) (6.2.6) 

and 

{Ur, Vr) (i/;, K) = ^~\'^r, K-) , (6.2.7) 

where is the inverse scattering operator deflned in Section 2.2 (in the case Uj. - Vr < 0, 
just replace X~^ with the identity). The flrst is a simple translation by the vector ?7fc^(t^) = 
Vkritr-i) J which is a function of Ur-i,Vj,r-i (see the Remark above). Therefore the result 
follows applying Lemma 2. ■ 
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7 Proof of the results 



According to the strategy of Lanford, once proven the uniform convergence of the two 
series (3.1.3) and (3.2.5) for short times, we shall conclude the validity results, namely the 
convergence of fj^{t) to fj{t), just proving the term by term convergence. Actually, by 
virtue of Proposition 2 in Section 5, it is enough to prove the term by term convergence of 
the series (3.1.4) to (3.2.5). 

In Section 6 wc have rephrased such expansions respectively in (6.8) and (6.2.1), i.e. 
sums over binary tree graphs of integrals over the (interacting or Boltzmann) backwards 
flows associated to the graph. Hence we must show convergence of the generic integral of 
this kind, 7^^(zj, t), to its analogue in the Boltzmann series, 7^^(zj, t). The present section 
is devoted to this problem. 

We stress once again the importance of the formulation of Grad (introduced in Sec- 
tion 3) which has been our starting point. In the language of Section 6 we could say that 
the terms in (3.1.3) that are absent in (3.1.4) collect all the interacting backwards flows in 
which two or more particles are created at some time ti (graphically, three or more lines 
emerge from a node of the tree). The use of reduced marginals has allowed to identify 
all these negligible terms and to isolate them from the contributions of order one, namely 
a^{j)T^^{zj, t). Now looking at (6.9), we see that these last object resembles very much the 
generic term in the series solution of the BBGKY hierarchy for hard spheres. Nevertheless, 
as we will explain in the next subsection, in the case of smooth interactions one has to be 
more careful in studying the behavior of 7J^(zj, t) for e small. 

7.1 The convergence problem: preliminary considerations 

Let us focus on 7^^(zj,t) and Ta-^i'^j-t). The integrand functions depend on the variables 
t„,i/„,Vj^„ completely through the trajectories of the IBF and the BBF respectively. In 
particular, the initial data f^j^n fo,j+n ai'e integrated over the time-zero configurations 
of the flows. Since foj^n converges to /oj+n by hypothesis, we must focus on the trajectories 
and prove that the IBF converges to the BBF for all values of t„, Vn, ^j,n outside a set giving 
a negligible contribution to the integrals. 

Looking carefully at the definition of C^("S) and Ci^) (see Sections 6.1 and 6.2), we 
realize that a great difference between them is generally caused by one of the following 
events: 

1. a particle (say j + i) created in the IBF interacts for a very long time (i.e. larger 
than 0(e)) with its progenitor; 
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2. a couple of particles {i, h) of the IBF undergoes a recollision, i.e. an interaction 
different from a creation (see the remark on page 29); 

3. a particle has a very large velocity, so that small differences between the two flows 
become large in a time of order 1. 

Item 1, which is obviously absent in the case of hard spheres, is controlled by cutting 
off the variables (z/j, f j+j) that lead to the singular scattering, and showing that they give a 
small contribution to the integrals. Here the main technical issue is an estimate of the time 
of interaction, such as that of Lemma 1 (or its generalizations in the cases of potentials 
with an attractive part discussed in Section 8). Similarly, item 3 is controlled by cutting 
off the energy of the system, i.e. the large values of |vj^„|. Item 2 is the most delicate. It 
requires to demonstrate that the contribution of recolliding trajectories is negligible in the 
limit e — > 0. 

To motivate our strategy in controlling the recollisions, we start by the heuristic analysis 
of one of the simplest non-trivial cases, namely that in Figure 4. At time t particles 1 and 




Figure 4: On the left, a simple case: r(j, n) = T{2, 1) = 1. The plus sign on the node recalls 
that 0"! = +. We want to estimate the contributions to the corresponding formula 7^(z2, t), 
coming from the recolliding trajectories of the IBF. An example of such a trajectory is 
symbolically represented in the figure on the right. 
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2 are in the final configuration Z2 = (a^^i, t"!, 3:^2, '^'2) £ M.2- We assume that the IBF is free 
up to time ii, when particle 3 appears with velocity at distance eui from particle 1, in 
outgoing ((Ti = +) coUision configuration. After the scattering between the couple (1,3), 
particle 1 collides with particle 2. This is a collision which is not a creation, i.e. what we 
called a recoUision. We shall imagine that Y — — Ci(^i) is order 1 while e is very 

small. We neglect the time of scattering between the pair (1,3) and approximate by Y the 
relative distance between particles 1 and 2 just before the scattering between 1 and 3. 

Denote by r/f the velocity of particle 1 between time ti and the time of the recollision. 
Then, the recollision implies a geometrical relation between W = V2 — f]'^ and Y. They 
must be chosen in such a way that there exists s G (0, ti) for which |.^|(s) — = ^- This 
is implied by the fact that W lies in the cone C(Y) with vertex 0. axis the direction of Y 
and tangent to the ball of center —Y and radius £, see Figure 5. Moreover, by the laws 

\ 

\ 

\ 

\ 








Figure 5: The recolUision-cone C{Y). 

of scattering (1.2) and (2.2.11), it is easy to see that 77^ belongs to the spherical surface 
of the ball centered in ^^li^ of diameter where V — — v^. In fact, fixed V and Vi 
(hence at fixed total momentum), 77^ moves over that sphere essentially as the scattering 
vector u (see Figure 1). In conclusion, jyf must belong to the intersection A of the cone 
V2 — C(Y) and the spherical surface described above. Clearly, at a given \V\, the surface 
measure of A is 0{e^) once assumed Y — 0{1) and e « 1. 
Now we want to estimate 

/ dvs du^u^- {v3-vi)l^-^^y , (7.1.1) 

J\V3\<R Jl^l-{V3-Vl)>0 

where the cutoff on ^3 has been added here to obtain an integral over a compact set. By the 
above discussion, it follows that a rather natural way to proceed is to express the integral 
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in terms of an integration with respect to V and u, so that we get 

J dV j^\\V\\V\'' j duB{u,V)l^^-^j,-^ (7.1.2) 

where V is the versor of V, we assumed also \vi\ < R, and B is the function resulting from 
the change of variables vi uj{vi, V) (see Eq. (1.3)). HE were bounded (as in the case of 
hard spheres) we would easily conclude that such a contribution is O(e^). Unfortunately, 
this is not true in many physically interesting cases, since B could not exist as a single- 
valued function and, even in each monotonicity branch of the scattering map p — > 0(p), 
it could diverge when the map becomes flat (see the Appendix). Thus to control the 
integral (7.1.2), we need to know properties of the scattering map (presence and strength 
of singularities), depending on the details of the potential. In this way it seems difficult 
to establish a unifled analysis of a large class of interactions (see the discussion in the 
Appendix). 

We propose a difi'erent method avoiding the use of the scattering cross-section (i.e. of 
the function B). This is based on two main ideas: 

- We work as much as possible on the Boltzmann flow, rather than on the interacting 
flow. Of course the BBF is much simpler since the interactions (= creations) are instanta- 
neous. Moreover, by virtue of the property described by the Remark on page 30, various 
parametrizations of the BBF, different from the usual in terms of t„, can be conve- 
niently used. In particular, the trajectories of the BBF can be parametrized by incoming 
collision variables. For these reasons, we find convenient to estimate the events in which 
some couple of particles of the BBF get closer than a certain distance (say on a scale 
shghtly larger than e). Indeed in the complement of this set the Boltzmann trajectories 
are close to the particle trajectories, as soon as the scattering time is small and the energy 
is not too large (which will be assured by an additional cutoff). 

- To estimate the above set of events, we use as much as possible the integration over 
time variables. From Figure 4 one can guess that, in general, only a small {0{e)) interval 
of values of will be compatible with the recoUision condition. 

However time integrations may produce singularities for special configurations of rela- 
tive velocities (see the Remark on page 41). Exploiting the global structure of the BBF, 
we will prove that such configurations are either excluded by the condition on the "initial 
datum" Zj e Qj, or they correspond to small set of values of relative velocities of incoming 
collisions, which will be estimated using the map of Section 6.2. 
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7.2 Proof of Theorem 1 



By the result in Proposition 2 of Section 5 and the reformulations of Section 6, the proof 
of Theorem 1 reduces to the proof of convergence of the generic term of the expansion, i.e. 



Proposition 3 Under Hypotheses 1-4, for all r(j, n),<T„ and {zj,t) e f2j x M+, 



(7.2.1) 



The aim is to apply the dominated convergence theorem to show that the trajectories 
of the IBF converge almost everywhere to those of the BBF. As already mentioned, we first 
need to "cut off away" pieces of phase space which correspond to trajectories of the IBF 
exhibiting recollisions, large scattering times, or high energies, and prove that they give a 
negligible contribution in the limit. Outside this properly defined set of "bad events", we 
will be able to estimate explicitly the distance between the interacting and the Boltzmann 
trajectories. 

In all this section and in the following we will keep fixed zj e [2j and t > 0. Moreover 
the times t„ will always be supposed to be ordered (see (6.1.1)), and the Un to satisfy the 
constraint implied by cr„ (Eq. (6.2.3)). In the present section we also fix T{j,n) and (T„. 

We start by focusing on the BBF <^(s) giving a new definition. Consider particle i 
and look at the graph of r{j,n). A polygonal path Vi is uniquely defined if we walk on 
the tree by going forward in time, starting from the time-zero endpoint of line i and 
going up to the root-point at time t. See for instance Figure 6. To Vi we may naturally 
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Figure 6: The fine closest to the dashed fine is the path Vi in the tree r(2, 6), with i = 8. 
The states of the particle associated to it via the BBF form the "virtual trajectory". 

associate a one-particle piecewise-free trajectory, built up with pieces of trajectories of 
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(different) particles of the BBF. More precisely, fixed a BBF with parameters (t„, i^n, Vj,„), 
denote ti^, - ■ ■ ,ti^, the (decreasing) subsequence of ii, • • • of the times corresponding 
to the nodes met by following the path Vi (rij being the number of such nodes, with the 
convention io = 0,tig — t): see the figure. We call virtual trajectory associated to particle 
i in the BBF, and indicate it by C{s) — (^'(s), 77'(s)) e with s e [0,t], the one-particle 
trajectory given by: 

{Ci(s) for s e \0,ti ) 
Cki,{s) for s e [ti^, ti^_J, < r < rij 

Observe that, during the time of existence of particle i in the BBF, C*(s) = Cii^)- 

Now consider a couple of particles (i, h) and compare their virtual trajectories. Call- 
ing "root" of Vi the root line of the tree to which Vi belongs, we have two possibilities: 
either the roots of Vi and Vh coincide (i.e. i and h belong to the same single tree), or not. 
In the first case, there exists (uniquely) a node of the tree where Vi and Vh merge. For 
any given couple {i, h) we introduce the subsequence of ti, • • • , t„ : 

t > t° > > > • • • > t"'" > = , (7.2.3) 

defined as follows. Time t^ is equal to t if Vi and Vh have different roots; otherwise, it is 
equal to the time (strictly smaller than t) of the node where Vi and Vh merge. The sequence 
^1 . . . ^ friih jg given by the ordered union of the times ti^, - ■ ■ , ti^. and th^, - ' ' ? th„f^ that 
belong to the interval (0, i°). Here riih is the number of such times, and has been put 

equal to zero by convention. See also Figure 7 below. 

We are ready to define a part of the "bad set" to be cutoffed. Let he 6 > 0. The set of 
S- overlaps, Af{S) C xR" x S^" x R^", is 

^f{S) = {tn,Jyn,^J,n min min If (s) - ^^s) | < ^ , (7-2.4) 

where t^ depends on the couple {i, h) under consideration. Notice that this set is completely 
defined via the BBF. Clearly, it depends also on Zj, t. 

Note also that the set J^{5) detects the 5-overlaps (namely when |f (s) — ^'^{s)\ < 6) 
of the virtual paths Vi and Vh excluding the time interval {t^, t°]. 

In the following, 5 > £ is a function of e going to zero as £ — > 0. Then, the first step in 
the proof is to show that the restriction of the integrals contained in T^^{zj,t) to the set 
Af{S) is arbitrarily small with e. To do so, consider the set of point-overlaps Af = A/'(0), 

M={tn,Un:y^j,n Hiin min [^(s) - ^^s) | = o| . (7.2.5) 

i<h se[0,ti] J 
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Obviously it is 

lim Ij^^s) = Iat ■ (7.2.6) 

We will now show that the set J\f has dA— measure zero. Precisely, we will show that the 
condition in (7.2.5) implies a certain number of relations between the integration variables 
that can be satisfied at most for a dA— null set of values. 
If we are in Af, then for some couple (i, h) there exists 

t* = max{s e [0, t^] s.t. \C{s) - t{s)\ = 0} . (7.2.7) 

It will be t* e [i'+^,t') for some I e {0, • • • ,nih}. I is the total number of interactions in 
the virtual trajectories of i and h between the overlapping time and For q — 0, - ■ ■ ,1, 
we define: 

Vl^v'is), for 5e(t«+\t'?), 

W = < - Vl ■ (7.2.8) 

We indicate by / e {0, 1, • • • ,n} the index such that 

f = tf. (7.2.9) 

Notice that cither f = t {f = 0) and ^ (because Zj G Qj) or t° < t (/ > 0) and 
y° = (because t° is the time of the node where Vi and Vh merge). A possible event in 
A/'(5) is pictured in Figure 7. 

Given a point in J\f, we observe preliminarily that we may assume: 

(i) if = 0, then ^ 0; 

(ii) / > 1; 

(iii) ^ 0. 

Assumption (i) corresponds to exclude subsets of the integration domain of codimension 3. 
In fact, by the laws of the two-body scattering, if F° = then ^ except for a single 
value of the velocity of the particle created at time i° = tf, namely vj+f must be equal to 
ijkfif^). Note that I — and y° 7^ is impossible because Zj e ilj. On the other hand if 
I — and y° = 0, then necessarily — 0, which we excluded by (i). Finally if VF' = 0, 
then the overlap takes place in the interval [i', t'~^) and this contradicts the definition of I. 

The point-recoUision condition is verified only if min^ \ Y^ — W^s\ — 0, which in turn 
imphes AW'' — 0, where W'' = Since 

y/ ^ yO _ ^ 1^9(^9 _ ^9+1) , (7.2.10) 
q=0 
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Figure 7: A symbolical drawing of virtual trajectories of two particles i, h in the BBF 
showing a 5— overlap after / = 5 interactions. Relative distances and velocities are indicated 
as defined in (7.2.8). In case A the two particles belong to different trees, while in case B 
they belong to the same tree, is the first (backwards) time when the particles reach the 
distance 5. 
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we have 

I 

= (yo - WH°) AW^ - Y^iiW^ - W"-^) A W^)]t'J . (7.2.11) 

q=l 

But, by the Remark on page 30, all the vectors involved in this relation do not depend on 
time. Hence as soon as [(1^^ — W'^~^) A W'')] ^ for some g, there exists at most one value 
of the time i« fulfilling condition (7.2.11). 
Otherwise, it will be 

W^^W'i = ^ for all g = 0,1,--- ,Z , (7.2.12) 

i.e. all (not vanishing) relative velocities are coUinear. In particular, Eq. (7.2.11) implies 

yo A = . (7.2.13) 

As said above we have two cases, which we treat separately. 

• Case y V 0, t° = 

Both y° and are coUinear with . Therefore 

yo A = , (7.2.14) 

which is excluded since e 

• Case yo = 0, t° < t. 

If all relative velocities are coincident, then no point-recoUision is possible since 
^ 0. Finally assume that U -.^W^ - W^'^ ^ for some ^ e {1, • • • , Z} and put 
[/ = . We have C/ A W^' = 0, which, together with AW^ ^ 0, implies 

UAW^ = Q. (7.2.15) 

Let — tfi, where f > f (recall (7.2.9)). We change variables 

(i^/, f^j+Z, ^f',mf') i^'f^ Vf, I///, Vf) (7.2.16) 

according to 
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In the first equation we used the map on page 30, while the second is a simple 
translation. We have Wq — and U — ±[77fc^,(i^,) — 77fe^,(ij/)] or U — ±[77^^, (i^,) — 
r]j+fi{tj,)]; see Figure 8. Prom the rules of scattering it follows that the vector U 
depends only on (i///, V//)^. Then Eq. (7.2.15) defines a subset of codimension two 
in the space of variables (V^, i///, V/'). By Lemma 4 of Section 6.2, the change of 
variables (7.2.16)-(7.2.17) is a one-to-one measure preserving map. Therefore the 
subset defined by (7.2.15) has dA— measure zero. 




Figure 8: Detail of the virtual trajectories of Figure 7-Case B, for times close to = tf. 
In the example af = +. The difference of relative velocities is U = — W^~^ = r]^ — r]~ 
(here /' belongs to Vi] if /' belongs to Vh then U — r]^ — ri~^), where rj'^ = rjk^X^f) 
r]~ can be equal to r]j^f'{tj,) or rik^.,(tji), depending on the structure of Vi- The variables 
describing the scattering are (z///, Wj+//), or alternatively {u'jr,, Vj,). 



So far we have proven that A/" is a null set and then it follows that the integral of 
any finite measure restricted to A/'(5) goes to zero with 6. But, in our Hypothesis 3, 
dA (Yl B^) foj+n is uniformly bounded by a finite measure for e small. Indeed, using 
conservation of energy, we can estimate it by 



dk 2 




e 2 



' ll/o,j+nll/3 1 



(7.2.18) 



where we used that the energy of the IBF at time t is purely kinetic if e is small enough 
(having fixed Xj outside the diagonals). Hence we have what we asserted before Eq. (7.2.5), 
that is 

/n 
dA{K, un, v,-,„) n B'i^i-, Vj+i - vUtmm)fSj+n{Cm = (7.2.19) 
i=i 



^It can he U = ±u){u; ■ Vf) ov U = ±[1//^ — (jo{uj ■ Vf)] where lj = V//) is the scattering vector. 
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for all Zj e Qj and alH > 0. 

Remark. In the above proof we try to use time variables, when possible, to show 
that the overlaps are rare. Otherwise we have to analyze geometrical conditions involving 
relative velocities of the virtual trajectories, namely Equation (7.2.12), which is proven to 
be a vanishing measure condition either by integrating in the incoming relative velocities 
of a node of the tree (case — 0), or by showing that the condition is impossible by 
definition of Qj (case y° 7^ 0). 

In the following section, when we will deal with quantitative estimates of the recoUision 
(or overlap) events, we will follow the same strategy, that is we will integrate in time unless 
we face the condition A W'^\ small for all g". Note that, in particular, this is the case 
of a sequence of central and grazing collisions in the virtual trajectories, for which the 
relative velocities may remain unchanged (remind that the scattering cross-section may 
possibly have concentrations on such collisions: see the Appendix). On the other hand, in 
this case the virtual trajectories are analogous to a free flow for which the recoUisions can 
be easily controlled. 

Besides H{5), we still have to take care of some additional subsets of the integration 
region. Let ji e (0, 1). Putting 

^1 " ^{f E£r^'<iiog6i} ' 

n 

= n ^{\i^3+r-Vl^{tr))f\yr\>et^} , (7.2.20) 

r=l 

a simple estimate as the one in (7.2.18) is sufficient to show that 

/n 
dA{tn,l^n,^j,n)YlB^{iyf,Vj+,-7ll{U)){l - 1^(5))(1 - lUI)/(^+n(C(0)) = . 
i=l 

(7.2.21) 

Thus, to obtain the final result we are left with the proof of 

/n 
1=1 

(7.2.22) 

Notice that up to now we did not use any property of the interacting flow but the conser- 
vation of energy. Now we have to examine in more detail the structure of the IBF and to 
compare it with the Boltzmann flow. Since we work in the complement of M{5) and in 
the sets in (7.2.20), we are actually in a favorable situation in proving that the distance 



41 



between the two flows is small. Indeed, as we will show in the following lemma, the IBF 
has no recoUisions and its differences with the BBF are only due to the scattering time, 
which is absent in the Boltzmann flow, and to the £— delocalization of the created particles 
(also absent in the BBF). 
Choose 

5 = £'-'*(log£)^ (7.2.23) 

Then we have: 

Lemma 5 // (t„,i/„,Vj„) is outside N^S) and inside the sets (7.2.20), then for e suffi- 
ciently small and all i = 1, ■ ■ ■ ,j + n 

max Us) - er(s)| < De'-^\ log£|i (7.2.24) 

se[o,t] 

for some D > 0. In particular, the IBF does not admit recoUisions. Moreover, r]%^{tr) — 
Vkritt) for allr ^l,--- ,n and 77? (0) = 77^(0). 

Proof of Lemma 5. The proof is based on a simple continuity argument. We proceed 
by induction on r proving that, for some D' > 0, 

lei(s) -er(s)l < D're'-''\\oge\'^ , s e (Wi,^r) (7.2.25) 

for z = 1, ■ ■ • , j + r, from which (7.2.24) follows since n is smaller then a constant times 
I logs I for e small. A byproduct will be that if particle i < j + r — 1, at time tr, has already 
completed the (possible) scattering with its progenitor (or its last son) in the IBF, then 
necessarily r/f (t^) = 7]i{tr) (from which the last assertion of the lemma follows, because the 
parameters are taken outside jV((5)). 

For r — the statement is trivial. Indeed, since we are outside A/'(5) with S > e and 
the states at time t of the BBF and the IBF are the same then, for s e {ti,t), (i{s) — 0(s). 

Let now s e {tr+i,tr) for r > and i — 1, • • • ,j + r. We consider the case (being 
the other cases easier) in which particle i in the IBF is interacting at time tr with another 
particle h. For e small and by inductive hypothesis, this can be true only if the two particles 
coincide with the couple {kr,j-\-r), or if they are a couple progenitor-son, {kr',j-\-r') with 
r' < r, that has not flnished the two-body scattering. Also, note that no other particle 
can be at distance < S/2 from i and h at time tr (for e small and n — 0{\ loge\)). 

Denote by s* e (tr+i,^r) the last (flrst backwards) recoUision time of particle i or /i in 
the IBF, that is |{^/(s) — ^f(s)| > £ for all h' ^ h,i and all s e {s*,tr), and same equation 
with i replaced by h. Then, for the same values of s, particles i and h behave as they 
were isolated and we have that |^|(s) — is equal to the same quantity evaluated in 
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max(s, tr — t), where t is the time of the two-body scattering (in fact, once the backward 
scattering finishes, it must be r]f{s) — r)i{s)). Hence, taking into account Lemma 1 and 
(7.2.20), we have 

lef (s) - m\ < IC^iU) - Utr)\ + 2^2/P\\oge\t 

< D'{r - l)e^-''\ log£|i + 2V2//3|log£|^£^-'' , (7.2.26) 



which implies (7.2.25) with D' = 2y/2j/3A, for s G (s*,t^). The same holds of course for 
particle h. 

In particular (taking e small so that n is 0(| log£|)) we have that |0(s) — 6 ("5) I ^ f up 
to the first recoUision time. Since we are outside A/'((5), for all h! ^ h,i 

\ew{s) - m\ > i^H'is) - m\ - m^) - e.'(«)i - ms) - m\ 

4 4 2 ^ ^ 

up to the first recollision time of i, h, h', and the same equation holds for particle h. But 
S/2 > e. Therefore, by continuity of the fiow, |C^,(s) - (s)| > | and \^l,{s) - ^l{s)\ > | 
for all times s e {tr+i,tr) and Eq. (7.2.25) holds in the full interval. ■ 

To conclude the proof of Proposition 3, we note that the above result, together with 
Hypothesis 4, can be also used to replace the initial datum f^j^^ ^y /oj+n in XTni'^j^^)^ 
that is to show that 

/n 
dA{K,u^,v,,n)YlB%u,;vj+, - rjliU)) |/i;.+JC(0)) - /o,,+„(C(0)) | = . 
i=l 

(7.2.28) 

Indeed we can restrict the above integrals for values of C^(0) in a compact subset of M.j+n: 
producing an arbitrarily small error. 

Using this last result, estimate (7.2.18) (with /oj+„ replaced by /oj+n), dominated 
convergence theorem and Lemma 5, Eq. (7.2.22) follows. ■ 



7.3 Proof of Theorem 2 

In this section we go through the steps in the proof of Theorem 1, estimating explicitly all 
the error terms arising in the limiting procedure, under the additional Hypothesis 5. We 
also need to take care of the dependence on n to guarantee the summability. For instance 
the bound (7.2.18) is useless in the present context because it grows as (n!)^/^. 
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In this section, for notational simplicity, we will denote by C all pure positive constants 
depending only on j3 and a. We also denote by £^2, • " " the errors in our procedure 
(i.e. it will be (z,,i) - /,(z,,i) = 

We start again noticing that it is enough to control the difference 

|f^(z„t)-/,(z„t)|, (7.3.1) 
since, by Proposition 2, the error due to this approximation is 

|^i| = |/f -/f|<C% (7.3.2) 

for t sufficiently small. 

Moreover, since 1 — Q!^(j) < jC^e^ iov n < N — j, we have also 

N-j 

n>o ro» <T„ 

A third simple error arises by neglecting the rest of the series expansion so that we 
focus on 

n 

Y.Y.ll\^L-r..\, (7.3.4) 

n=0 r(j,n) cr„ 

where 

n = D"\ \oge\ (7.3.5) 

with £ small enough. The error generated by this last cutoff is bounded by the remainder of 
the geometric series appearing in formula (5.10), therefore if we choose D" = \ log(tC^ q,)|~^ 
it is 

\^^\<Y.T.T.i^^r. + r..)<C^e, (7.3.6) 

n>n r(j,n) <Tn 

for t sufficiently small. 

Before considering the difference (7.3.4), we estimate the error caused by the restriction 
of the integrals to suitable sets of integration variables, where the convergence of the flows 
(Lemma 5) is guaranteed. 

However, flrst of all, to simplify the expression of the integrands avoiding estimate 
(7.2.18), we get rid of (H B^) performing the decomposition {(H B^) > e~^} and {(H B^) < 
e~^} (and the same for (fl B)) for a suitable A e (0, 1). Moreover we shall often use in the 
following the bounds 



/oj+n(C(0)) < C^+"e-^n- , (7.3.7) 
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consequences of Hypotheses 2 and 3, the conservation of energy and the fact that the 
energy at time t is purely kinetic if is outside the diagonals and e is small enough. 
We need to estimate 

n „ n 

n=OT{j,n) (Tn «=1 
P n 

+ J o?A(t„, Vj-„)l|(nij)>£-A} Yl B{ui; Vj+i - %,(t+))/oj+„(C(0)) 



< £^ E E E / [ (n ) + 



(n 



B 



(7.3.8) 



Remind now the expression of and Eq. (6.4). Since X^fc^i^ {vlMi))'^ bounded by 
the total energy Yli=i '"1 (being the potential positive) , it follows easily that 



j+n 



E(ns') <2"n o+^K+i+E"? 

r(j,n) i=l \ 1=1 

The same estimate holds for X^rQ n) (11 -^)^ ■ Therefore 

n>0 i=l ^ 



(7.3.9) 



where we used the bound 



i+n 

E 

i=l 



1 --e-T-'+"7;2 ^ 



.i+i„2 ^ 4n 
e5 



The integral on the velocities in (7.3.10) factorizes so that 

i^4i <e'E^'^"^(^'+^)"- 



n>0 



Since 



(j + nY ^ {j + n) 



ni {j + n)\ 

we have that (7.3.12) is bounded by a geometric series. Hence we conclude that 

1^4 1 < C^e^ 



(7.3.10) 



(7.3.11) 



(7.3.12) 



(7.3.13) 



(7.3.14) 



for t sufficiently small. 
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At this point we just follow the hnes of the proof of Proposition 3. Using the notations 
introduced on pages 35-36 to define the set of 5— overlaps in the BBF, we want to estimate 



cN 



i^si < E E E [ / ^^ i{(n^^)<e-} (n^O 

+ JdA l{(UB)<e-^} (n^) Ms)foj+n] • (7.3.15) 

By (7.3.7) it is 

l^^l ^ ^"'E E E^'^" / ^MiS)^'^"^ ■ (7-3.16) 

n=Or(j,n) CT„ 

We have the following crucial result: 

Lemma 6 Let 5 = £^~'*(log£)^ (see Eq. (7.2.23)j. Given zj e if e is sufficiently small 
and 1 < n <n, 

J rfA(t,,i/,,v,,,)l^(5)e-inii"-l < C^-+"^5i|log£|i . (7.3.17) 

Here the choice of 5 is not strict and is determined by Lemma 5. 

Proof of Lemma 6. First notice that, given a point in N{5), there exists 

tl = max{s e [0, t^] s.t. |f (s) - i^{s) \ < 5} (7.3.18) 

for some couple of particles (i, /i). Substituting t* by t|, we may define /, F^, W^, f as after 
(7.2.7); see also Figure 7. 

The case / = and = 7^ is made impossible by Zj G as soon as 5 (i.e. e) 
is smaller than a constant depending on Zj. Conversely, the case / = and t*^ < t, y° = 
occurs whenever a creation in the BBF is such that the two particles progenitor-son do not 
separate enough before their next (backwards) interaction. In formulas, |iy°|(t° — t^) < 5. 
This case is controlled by introducing 



- H ^{\{vj+r-'nUr{ttmtr-tr+^)>2U'^/^-} (7.3.19) 

r=l 

(remember that the modulus of relative velocity is conserved at collisions) . Clearly it is 

n 

^M{6) < 1o1a/-(5) + E ^m+r-vkAt+mtr-tr+i)<2tsv^^} ' (7.3.20) 
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■Jj+r — 'lkry 

r=l 



The reason for the choice of the threshold in definition (7.3.19) will be clear soon. For the 
moment note that we have 

J2 J ^A(t„,i.„,v,,„)l{|(,^^^_^^^(,+))l(,^_,^^^)<,,,./5}e-^S^-"^? < C^^^—^S'^ . (7.3.21) 

In fact, performing first the integration in dtrdvj+r, reminding that r]k^{t'^) — r]k^{t~_i) is 
independent on tr, Vj+r, setting s — tr — tr+i and V — Vj+r — Vkri'tt)^ 



tr-l 



/3 2 



tr+1 

< ds dvj+re-^''^+r + ds I d\V\in\V\'^ 
Jo J Js^/^ Jo 

< C5i . (7.3.22) 

Since the integrals in the other variables give C"~^i"~^/n!, we get (7.3.21). 

Suppose now that the considered point in N{5) is such that l^lj^f^s) — 1. The 5— overlap 
condition is verified only if 

min - W^s\ < 5 (7.3.23) 

s 

for some s G [0,t'), with given by (7.2.10). It must be / > 1 and ^ 0. Moreover, 
the relative velocities cannot be all too close to each other, i.e. the two characteristic 
functions imply 

I 

^\W'i -W'^-^\> si . (7.3.24) 

9=1 

Otherwise it would be \ W'^ — W'^\ < 5^ for all g, thus using (7.2.10) and (7.3.23) we would 
deduce that 

\Y^ -W^s\<5 + 5h (7.3.25) 

for some s > (i° — i'), which is forbidden, for 5 small, either by Zj e Qj (case Y^ ^ 0) or 
by definition of (case y° = 0). 

Eq. (7.3.23) imphes \Y^ A < S, where = Using again (7.2.10) we have 



5 > 



q=0 



{Y° -Wh°) AW^ -Y^liW^ -W^-^) AW^)]t'^ . (7.3.26) 
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Since the vectors involved in this relation do not depend on the times (Remark on page 
30), we can exploit the integration in the variables to estimate the set defined by this 
condition. However, a singularity will arise when the vector in the square brackets is small 
for all q. Let us focus first in this case, which is the most delicate. Assume that 

I 

J^liW^ -W^-^) AW^\ <5^ . (7.3.27) 

q=l 

Notice that if Eq. (7.3.27) is not satisfied, then \(Wi* - W^*-^) AW^\ > 6l /I for some q*. 

Again we may infer that condition (7.3.27) is not possible in the case = t, F° 7^ 0. 
Indeed, the above inequality trivially implies 

\W^ AW'\<5^ , (7.3.28) 

therefore (7.3.26) gives 

\Y^ AW^\<6 + 2tS^ . (7.3.29) 
Putting together the two last equations we have 

\Y^ AW^\<C {\Y^\ + \W°\)SI , (7.3.30) 

which is excluded, for 5 small, by Zj e f2j. 

Let us study (7.3.27) in the case t° < y° = 0. By (7.3.24), there exists a ^ e {1, • • • ,1} 
such that 

= - W^-^ (7.3.31) 

has modulus 

\U\> — . (7.3.32) 

t 

But \U AW^\ < by (7.3.27). Then putting [/ = we have 

\UAW^\<n6'^. (7.3.33) 
On the other hand, by (7.3.28), either 

< 5i , (7.3.34) 

or A W''\ < 6^ which, together with (7.3.33) and the constraint n < n = 0{\ log£|), 
finally gives 

AU\ <C6'^\ogs\ . (7.3.35) 
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We will use this formula to estimate the considered events, taking advantage from the fact 
that U depends only on the impact vector and relative velocity describing the interaction 
occurring at time in the BBF (as already pointed out on page 40 and Figure 8), and 
that is the (incoming) relative velocity of the interaction at time t^. 

We shall summarize the discussion above as follows. Denote K and respectively the 
outgoing and incoming relative velocities of the collision at time tr in the BBF. If — tr, 
we use the notation Uj. — U^. This is a function of v^, K only. We have 

riif^ I n n n 

itl^iS) ^J2J2Y1 ^^!f{S) + h\Vr\<sy'} + ^{\v}AUf,\<cs^/^\iogs\} (7.3.36) 

i,h 1=1 q*=l ' r=l f=l /'=/+! 

where l<i<h<j + n and 

i^n,^j,n the virtual trajectories of i and h satisfy (7.3.26), 

with |(iy«* - W'-') A > 5l/l^ . (7.3.37) 

Once fixed all the variables but P* , if Equation (7.3.26) is verified, then t** belongs to 
an interval of length smaller than SliW^* - W*''^) A W^\-\ If we are in X-f{S) this is 

2 

bounded by nSs, so that 

J dA(K, u^, v,-„)l^,.*(,^e-i S^=>? < C^-^^^sl . (7.3.38) 
Changing variable Vj+r — >■ K we easily find 

J ciA(t„i/„v,-„)l{|v-^|<52/5}e-i^'i"-' < C"^5f . (7.3.39) 



Furthermore, it is 



/ 



dK{tn.Un,^j,n)l{\v>.U,,\<C5y^\loge\}^ i "I < C"^5i|log£|i . (7.3.40) 



To prove this last inequality, it is convenient first to introduce a further restriction to 
the set {f |V/'p < 4| log^l} where Vfi = (fj+/' — (^//)). If the opposite inequality holds, 
then either 1^;^+/' | or \rjkj,, {tf')) \ cannot be smaller than \/2//3\ log£:|, hence the total energy 
must be larger than l/^|log£|. Therefore, using the energy-cutoff If defined in (7.2.20), 
the error produced is bounded by 

dA(t„,i/„,v,-„)e'^^-"^? < C"-£5 , (7.3.41) 
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which is in turn certainly smaller than {Ct)'^/n\55, being 5 given by (7.2.23). We are left 
with 

J dA{tn,U^,Vj,n)^y,^u^^^^cS^/,^,^s^^^^ ■ (7.3.42) 

We change the integration variables according to (7.2.16)-(7.2.17). Note that 



e-l^hfe-^^^f = e-K''/-'"^("^-''/)+^^/^*/0'e-Kv+'?V(*/'))' , (7.3.43) 



where o^/ = a;(i/j, V^) is the scattering vector at the collision. Since Vf varies in a compact 
set, we bound the second exponential simply by 1, while the first exponential is estimated by 
g-2 ^ ^jigpg r]f;^(t+'^ depends only on the variables with index strictly smaller 

than /. Performing first the integrations in dVfdVf we find 

( ^^^'^ \ 

<C 1+ ^ (5i|log£|i . (7.3.44) 

Integrating in the remaining variables we readily get Equation (7.3.40). 
Collecting all the estimates, Lemma 6 is proved. ■ 

Substituting Eq. (7.3.17) in Eq. (7.3.16), performing the sums and using (7.3.13), we 
conclude that 

l^^sl < C%-^5t|log£|i < el--^^'^\\ogEp (7.3.45) 

for t small enough. 

We turn now to the estimates of the errors coming from the truncations defined in 
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(7.2.20). Proceeding as above we have 

i^ei ^ E E E [ / ^^ i{(n^^)<-n (n^^) (1 - (1 - nil) 

n=oro>) o-„ 

n „ 

< C%i-^ + E E E dA {1- 1|) e-^ ^'i"''' , (7.3.46) 



n=0 r(j,n) CT 

where in the last step we used (7.3.41). Moreover 



[dA{l- 11) e-i <Y.ldA l(|(.,,.-.^^(..))A..|<e.}e-' ■ (7-3.47) 

r=l 

We perform first the integrations in dvj+rdi'r- Setting = {''^j+r ~ Vlri^r)), the angle 
between K and Ur and using the parametrization Ur — > (p, (p) where p — sin a and e 
[0, 27r), we have 

//3 2 
dVj+rdUrl{\(v^^^-r,l^(tr))Au,\<e>^}e~^''^+'~ 

e** 1 — 

< / civj+^(27r)2 / dp ^ e~2^^- + (27r)2 / dp—^= [ " d|K|47r|K 

J 7o V 1 - p2 JeM >/l - p2 Jo 

< Cs^'^ . (7.3.48) 

The integrals in the other variables give C"~^i"/n!, therefore performing the sums as above 
we obtain 

l^el < (e^"^ + ■ (7.3.49) 

Now we shall estimate what is left of the differences in (7.3.4). This gives two errors: 
one is due to the convergence of the initial data (formula (4.7) in Hypothesis 5) and the 
other is due to the convergence of the IBF to the BBF (Lemma 5). The first is 



2 



n „ 



-"-2 



n=0 r{j,n) <Tn 



foj+niCm-fo,j+n{Cm . (7.3.50) 
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Since we integrate outside Af{S), the BBF satisfies C(0) G Mj{S). But 5 = £^~''(log£)^. 
Thus, applying Lemma 5, the IBF must satisfy C^(0) £ •M.j{e) for e small enough. Hence 
Hypothesis 5 together with conservation of energy lead to 



n p 



7 A _^T^+",,2 



"=0 r(j,n) o-„ 

having performed the sums in the usual way. 
Finally, the last error is 

n „ 
n=OT{j,n) cr„ 



(7.3.51) 



/o,+n(C(0))-/o,,+n(C(0)) 



-"-2 



(7.3.52) 



Lemma 5, the regularity assumption (4.8) in Hypothesis 5 and conservation of energy at 
collisions imply that we can bound the modulus in the last line by 



< C^'+"e- 2 Efir ( j + n) ^ L>£^-'^ I log £ I i , 

(7.3.53) 



for s sufficiently small. Therefore, proceeding as above we have 

l^sl < CV-'^-^llogeji . 



(7.3.54) 



Putting together all the errors £i, - ■ ■ ,Ss and optimizing on /j, and A we conclude that 

1 



cN 



(zj , t) - fj (zj ,t)\<C^e^ for any 7 < 



6 



(7.3.55) 



8 Stable potentials 

In this section we show how the techniques used in proving Theorem 1 can be extended to 
treat a fairly larger class of potentials, including those with an attractive part. 

The potentials $ considered in the present section satisfy the following conditions. 

Hypothesis 1' $ = ^{q),q € is radial, with support \q\ < 1. We further assume 



52 



1) either $ G C^iM?), or $ G ^^(M^ \ {0}) and $(g) ^ +00 as q ^ 0; 
$ is stable. 

In what follows we will use the usual notational inconsistency $(r) = $||g|=r. 

We remind (see e.g. [14]) that an interaction is stable if it fulfills the following condition: 

f/(gi, ... ,q^) = Y, <f (ki - qh\) > -Bj (8.1) 

i<h 

for some constant B > Q . In particular, $ is positive (possibly diverging) at the origin. 
We also remark that condition (8.1) ensures the existence of the Partition Function and 
hence the existence of an equilibrium measure. In our context the stability will be used to 
guarantee that Hypothesis 3 implies the bound f^^ < c^e'^^^""^ (where c = e"^^^), which 
is crucial in our proof. 

The potentials $ we are considering include a sort of truncated Lennard- Jones potential 
(see Fig. 9 below). 




Figure 9: A cutoffed Lennard- Jones potential. 

We note that the proof presented in Section 7.2 depends on $ only through the scat- 
tering time r=i,. Therefore, in trying to extend it to the present situation, the crucial point 
will be the control of r^,, in absence of Lemma 1 of Sec. 2.2 which is not valid anymore. 

According to Section 2.2, we reduce the two-body particle system to the central motion 
of a single particle, with mass | and velocity V (the relative velocity of the two interacting 
particles). 
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We recall the formula yielding the interaction time: 

n^V2 [ dr— — , (8.2) 

where |y| > is the modulus of the relative velocity before the collision, p is the impact 
parameter, L = \pV\ G [0, \V\] is the modulus of the angular momentum and r* is the 
infimum of the distance from the origin during the scattering process, r* is given by 

n = max ja; G [0, 1) s.t. Y ^ 2^ ^ ^^^^'^ J " ^^"^^ 

Before establishing the following lemma in which we control the scattering time, we 
discuss the new difficulties we face in presence of an attractive part. Consider for instance 
the potentials described in Fig. 9, with a single negative well. The effective potential 

*(^) + ^ - T ^^-^^ 

may have two critical points, and (minimum and maximum respectively), when L 
is sufficiently small; see Figure 10. Fixing a value of L for which such critical points do 



Figure 10: When the (real) interaction has attractive parts, the effective potential (drawn 
in figure) can have local maxima for a given value of L. 

exist, there are values of V for which 

— ^ 2$(rM) + (8.5) 

In this case the trajectory is close to an unstable periodic orbit and is very large. The 
two particles turn around each other many times and remain trapped for a long time. 
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Clearly such situations are pathological and must be excluded in order to have a kinetic 
picture. Actually the following lemma says that the set of such pathological events has a 
small measure, although we do not give explicit estimates. 

Lemma 7 Given rj e (0, 1) and K e M+, there exists a set B{rj) of pairs {u, V) such that 

dV\u ■ V\lB(n) ^ (8.6) 



[ du [ 



my) + ■ (8.9) 



•y<o} 

as 77 — >■ and such that, for {u, V) ^ B{ri) and \ V\ < K, 

n{u,V)<cp{ri,K), (8.7) 
where ip{r], K) is a positive function which may possibly diverge as rj ^ or K ^ +00. 

Proof. We can easily find the set of pairs (L^, y^) for which r* diverges. Such a set is 
included in the set for which there exist local maxima of the effective potential (8.4). The 
critical points satisfy 

^'{y) - ^ = . (8.8) 

Therefore the pairs (L^, V'^) corresponding to a divergence of the scattering time r* must 
satisfy (8.8) and 

Yl 

2 ' 2y^ 

This last condition is due to Eq. (8.3), while (8.8) ensures that the orbit reaches = y in 
an infinite time. 

Consider now the curve C in the plane, y e (0, 1) (X, F), whose parametric equations 

are 

' X = 2^'(y)y'^ 

^ ' . (8.10) 

Y ^A^{y) + 2^\y)y 

Then the set of singular values of (L^, V'^) lies inside the restriction C of this curve (X, Y) 
to the "physical" subset 

|(X,F) s.t. y>0, 0<X<Fj. (8.11) 

Clearly when $ is bounded the curve C is extended by continuity to y = (for which 
= X(0) = 0, = y(0) = 4$(0) are indeed singular points of r*). Note that, when $ 
is unbounded, the parameter y spanning C is bounded away from zero, since it cannot be 
smaller than ro := min{x e (0, 1] s.t. ^'{x) > 0}. 
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Denoting by B{{X, Y);r)) the disk of center {X, Y) and radius rj, we introduce the tube 

r{ri)^ U B{{X,Yy,rj) (8.12) 
(x,Y)eC 

and its restriction to the physical region 

f{v)= U B{{X,Y);n). (8.13) 

{x,Y)ec 

Now observe that, due to the smoothness of the set C has finite length so that 

\nv)\<Cv, (8.14) 

where \A\ denotes the Lebesgue measure of the set A. 
Consider the set 

f(v)^B{{0,Oy,rj)U{\V\>K}. (8.15) 

Its complement Q{ri,K) in the physical region (8.11) is relatively compact. Therefore, by 
continuity of t^{L'^, V^) in the set (8.11) deprived of C (and hence in the closure of Q), we 
have 

n<ip{ri,K) (8.16) 

in K), for a suitable positive function (/?, possibly diverging as — ?> or i^T — ?> +oo. 

To prove the required continuity of r*. we observe first that = r^lL"^, V"^) is continuous 
outside C. Fix a point (L^, V^) ^ C in the set (8.11). Then, for any 7 e (o, IzIiMjYol^ ^ 
the integral 



V2 / dr ^ (8.17) 

is continuous in (Lq, Vq), being the integrand bounded. On the other hand, 

V2 / dr ^ — > (8.18) 

as 7 — )■ 0, uniformly for {L'^,V^) G -B((Lq, Vg^); 5), if S is small enough. This can be seen 
by an argument as the one in Lemma 1, namely using the estimate (2.2.8) replacing the 
integration interval (r*, 1) with (r*,r* + 7). 

To conclude the proof of the lemma, we introduce the set 

B{r]) = {{u,V) eS'xR'\ (L^y2) e f{v)UB{{0,0);r])} (8.19) 
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where L — \v /\V\. Setting cos a — v - V^V — V/\V\ and noticing that L — \Vsm.a\, the 
left hand side of (8.6) is bounded by 

n n POO POO pit 

du dV / (i|F||V^|^| cosq;|1b(j,) = 87r^ / o?|V^||V^|^ / dasinal cosq;|1b(^) 
Js^ Js^ Jo Jo Jo 



< Att^ / dV^ / dLHe^n) < Crj 
Jo Jo 



(8.20) 



for T] sufficiently small. ■ 

We are now in a position to establish and prove the main result of the present section. 

Theorem 1' (Improved) Under the Hypotheses V and 2-4 of Section 4, there exists to > 
such that, for < t < to and j > 0, 

lim ff{t) = f,{t) (8.21) 

Ne^=l 

uniformly on compact sets in Qj. 

Proof. We just mention where the previous proof of Theorem 1 requires modifications 
and how to do them. 

The proof consists of two parts, namely the short time estimate and the term by term 
convergence. 

As regards the short time bound in Section 5, by virtue of the stability property, it is 
natural to modify the definition of the Hamiltonian by setting 

//s(z,.) = /7(z,.) +jB>\Y.v^,>Q. (8.22) 

i=l 

Consequently we introduce the norms (5.1) replacing H by Hb- Next we deduce estimate 
(5.3) by using (8.22) and the fact that Hb satisfies the inequahty (5.5), i.e. 



^ j+l+m 

HB{Zj+l+m) = Hb{Zj) + HBiZj,i+m) > Hb{Zj) + 3 E (^'^S) 



i=j+l 



We use the stability of the potential only in this part of the proof. 

Now we pass to analyze the term by term convergence. Everything is going on as in 
Section 7.2, with the only difference that we replace 



n 



^2 = n ^{\iVj+r-Vl{tr))AUr\>e^} (8.24) 



r=l 
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by 

n 

I2 = n ^{i''r,Vj + r-Vl^{tr))mv),\vj+r-Vl^{tr)\<K} ■ (8.25) 

r=l 

According to the form of the function (/?, we can choose rj — rj{e) and K — K{e) such that 
77 — > and K ^ 00 with e and the scattering time of the coUisions associated to the nodes 
of the tree (in macroscopic variables) is bounded by e(p{r], K) < Ae^~^ whenever 1| = 1 
(by virtue of Lemma 7). Therefore Lemma 5 (hence Eq. (7.2.22)) is still vahd and the 
proof is completed by observing that, by Lemma 7, ^(77) is a set of vanishing measure (so 
Eq. (7.2.21) holds). ■ 

9 Concluding remarks 

We conclude by discussing some additional remarks. 

1. The potentials we have considered are fairly general, but the basic hypothesis is the 
short-range assumption. 

Prom the very beginning of the Kinetic Theory, Boltzmann himself (see [1]), following 
Maxwell [12, 13], considered only inverse power law potentials, besides the hard-sphere 
system, originally investigated in deriving his famous equation. This is probably due to 
the good scahng properties of such potentials. Moreover the differential cross-section is 
well defined, even though the total cross-section is diverging because of the long range 
of the interaction. On the other hand, it is not clear whether the Boltzmann equation 
associated to these potentials can indeed be derived under the low-density scaling. 

A simpler problem would be to consider a sequence of potentials with range gently 
diverging with N. This problem eludes the present techniques so that we consider it as an 
interesting, open problem. 

For an analysis concerning the much easier problem of the validity of the linear Boltz- 
mann equation for long range potentials, see reference [4] . 

2. In the present paper we give an explicit estimate of the error in case of a completely 
repulsive potential (Theorem 2), while, for stable potentials, we only show the convergence. 
It would be interesting to develop a constructive proof of convergence also in this last case. 
This would require a more precise estimate of the scattering time to improve Lemma 7. 

3. The present validity results, as the ones in the previous literature, are formulated in 
a canonical context, namely, for any e > 0, the number of particles N is authomatically 
fixed. An equivalent formalism is the grand-canonical one. Here the number of particles 
is random but the density is fixed. 



58 



More precisely consider, for a given £, the phase space of the system as 

M=[jMN (9.1) 

Af>0 

where A^at is the A^-particle phase space (see (2.1.2)). For z ^ Aiwe define the dynamical 
flow by solving the Newton equations in each A^jv- Similarly we define a symmetric 
probability measure on A4 by means of a sequence of symmetric probability measures 
in each Mn' 

W'lM.^e-^^^W''. (9.2) 
The sequence {gj}^i is defined by 

,.N 

^I = E--"|r^f (9-3) 

N>j 

where g^^ are the marginals of . Therefore gj{zj) are the probability densities of finding 
the first j particles in Zj. Their normalizations are 

,.N 

N>j 

which is the probability of finding more than j particles. Then one defines the reduced 
marginals accordingly and it is easy to derive the equivalent of the Grad hierarchy for 
them. 

Note now that the average number of particles is (N) = ji^. Therefore the low-density 
limit will correspond to /Xg — > oo and e^^^ > 0. 

It is not difficult to realize that the validity result of this paper can be formulated and 
proven also in this context. 

4. We have considered the particle system in the whole space. If we want the system 
to be confined in a bounded box, we have to specify the boundary conditions. Assuming 
specular reflections, there are additional difficulties which we have to overcome. First the 
dynamical flow is only almost everywhere defined (see [11]), but this (as for the hard- 
sphere systems) does not create real difficulties. However the analysis of the recoUisions 
requires some extra geometrical arguments. 



59 



A Appendix (on the cross— section for the Boltzmann 
equation) 

In this appendix we give sufficient conditions on the interaction for having a single-valued 
differential cross-section (and we show some counterexample). We also study the bound- 
edness properties of the cross-section. The issue is relevant both to motivate our strategy 
and to know whether S is a well behaved single-valued function in the usual form of the 
Boltzmann equation, Eq. (1.3). 

The assumptions on $ are those estabhshed in Hypothesis 1', but possibly allowing a 
discontinuity of the first derivative at \q\ — 1. 

Consider the planar scattering process of a particle of unit mass. We use the notations 
of Section 2.2 and of Figure 1. In particular, we denote by p the impact parameter (by 
symmetry we may focus on < p < 1) while the scattering angle is x = tt — 2© and the 
energy in the laboratory > 0. 

The differential cross-section is defined through the map p — p{Q, \ V\), by 



2|sin(2e)| 



dp 



(A.l) 



Therefore we need to analyze invertibihty of the map 6(p). 
The classical formula for © is 



©(p) = arcsin p + p / dr , (A.2) 

where r* is the minimum distance of the central motion from the origin, satisfying 

1-^-^ = 0. (A.3) 

For purely repulsive potentials with a singularity at the origin, the limiting values are 
0(0) = and 0(1) = 7r/2. In general, it is = nvr + 0' for some 0' e (0,7r), where n 
is the total number of counterclockwise turns that the trajectory makes around the origin 
(see Section 8). 

While the first term in the right hand side of (A.2) is an increasing function of p. the 
second term is clearly non monotonic (in fact it goes smoothly from to when p — )> 
or p — 7> 1 and hence r* 1). Following [4], we set y = p/r and perform the change of 
variables 

2$(2) 

-^ + y2 = sinV, (A.4) 
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to get 

^/ \ • r^"^ , sin 09 / A r\ 

e(p) = arcsinp+ / d^p ^jTTzy • (A.5) 

^arcsinp y - 

The advantage of this formula is that the integrand is not singular in the integration region 
and we can easily compute the derivative with respect to p. 
A straightforward calculation leads to 

\ ( ^ I . 

1 - . .^n-. (A-6) 



,2 



dp .JY^^ \ 

f7r/2 



+ dip ^ 

Jarcsinp (y__|L_$/(£; 



V^y^ \y J V^y \y J V^y^ \ \y 



for < < 1, where indicates the limit of the derivative as \q\ — )■ 1 from below. 

In formula (A. 6) we are also considering the case in which $ has a discontinuity of the 
first derivative in |g| = 1 as it is the case of the inverse power law potential restricted to 
the unitary interval treated in [4]. However, for the case of smooth potentials as the ones 
considered in the present paper, the first term in the right hand side of Eq. (A. 6) is absent. 

The following considerations can be deduced from Eq. (A. 6). 

1) The ratio p/y ^ g{(fi) as p — )■ 0, where is a positive function of (f which form depends 
on $ and V'^. Then the extremal values of our derivative are: 

*'<^"'^° . (A.7) 
$'(1) = ^ ^ 

2) The monotonicity property ^ > translates in a quite complicated condition on the 
potential $. A convenient sufficient condition is given by the following assertion: In the 
considered class of potentials, if for all q with \q\ e (0, 1) 

\q\<^"{\q\) + 2^'{\q\)>0, (A.8) 

then ^ > for all p e (0, 1),V^ > 0. This condition is derived also in [15]. 

Condition (A.8) can be easily checked for a large subset of potentials. For instance 
any potential of the form $(g) = — 1^ S\q\<i, k > 1, satisfies the condition, hence has 
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strictly monotonic map. Cases which are smooth in |g| = 1 can be constructed from the 
previous by using a smooth junction. For instance^ 




(A.9) 



Figure 11: Map Q{p) for the potential given by Eq. (A.9), with 6 = 0.l,k = 20 and Eq = 
= 9. 

3) The monotonicity property ^ > is in general not true when condition (A. 8) is 
violated. A first example is any smooth, positive, decreasing and bounded potential (for 
which G(0) = 0(1) = 0, that implies the existence of at least two monotonicity branches). 
We give two different examples of potentials singular at the origin. 

• Formula (A. 6) indicates that the sign of the second derivative of $ is relevant when 
we ask about monotonicity of the map. In fact, examples of non monotonic maps 
can be constructed when $" is not always positive, for instance by taking $ very 

•^Observe that this is a function C^{R'^) with a jump in the second derivative for [gj = 1 — S. The 
parameters k and S can be arranged in order to ehminate this discontinuity (e.g. S = 1/10, fc = 71). 
Nevertheless, all our discussions are still valid if in the initial assumptions on the potential we require that 
$" is just piecewise continuous and bounded outside any ball centered in the origin. 
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close to the characteristic function of |g| < 1. If we consider the function 

= -etan ^^arctan- + 1^1 ~ + ^ ' (A. 10) 

numerical simulations show that the map Q{p) is non monotonic for e << 1 as shown 
in Fig. 12. 



E=0.1 



0.2 0.3 0.4 0.5 

P 



0.6 0.7 0.8 0.9 



Figure 12: Map e{p) for the potential given by Eq. (A.IO), with e = 0.1 and Eq = V"^/! = 3. 

• Even when $" is nonnegative, the mapping can be non monotonic if Eq. (A. 8) fails, 
an example being 

m = < 5(1 -\q\) 6<\q\<l, (A.ll) 

[ |g| > 1 

We checked numerically the non monotonicity of 0(p) in the case S = 0.1, k = 4 (see 
Fig. 13). Another example similar to the previous one but with continuous derivative 
in |g| = 1 can be constructed again by using a smooth junction. 

4) Any time the condition ^ > is violated, we have a singularity of the cross-section, 
even though single- valued. In particular, if $'(1) = (i.e. the force is smooth) and Q{p) 
is strictly monotonic, we still have a divergence of ^ for G near to 7r/2 (p = 1), that is 



sin(26)cr$ 



+ 00 



(A.12) 
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6=0.1 



qI 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

P 

Figure 13: Map @{p) for the potential given by Eq. (A. 11), with 6 = 0.1, k = 4 and Eq = 
= 0.5. 
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